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MACH l^EFLECTION, MACH DISC, AND THE ASSOCIATED NOZZLE FREE JET FLOVfS 

I-Shih Chang, Ph.D. 

Department of Mechanical and Industrial Engineering 
University of Illinois at Urbana-Champaign , 1973 

The numerical method involving both the method of integral relations 
and the method of characteristics has been applied to investigate the 
steady flow phenomena associated with the occurrence of Mach reflection 
and Mach disc ft'om nozzle flows. The solutions of triple-shock inter- 
section have been presented in detail. The regime where Mach configur- 
ation appears is defined for the inviscid analysis. The method of inte- 
gral relations originally developed for the blunt body problem is modi- 
fied and extended to the attached shock wave and to the internal nozzle 
flovj’ problems. This method, together with the method of characteristics, 
is- -adopted to study Mach reflection from two-dim.ensional -overexpanded noz- 
zle flows Vfhereas the axisymmetric underexpanded nozzle flows with Mach 
disc are studied by using the method of characteristics exclusively. The 
reflected shock configuration, as v?ell as the accompanying flow field in- 
cluding the Mach stem height in two-dimensional cases and the Mach disc 
radius in axisymmetric jet flow, is established through the consideration 
that the central core flow downstream of the Mach -configuration should 
reach a state which is equivalent to choking for a maiform one -dimensional 
flow. In the axisymmetric underexpanded nozzle study, it is found that 
the rotationality x'esulting from the entropy difference must be considered. 
The numerical results are. compared vrith the experim.ental data of Love and 
Grigsby for the axisymmetric gas jet. An approximate method based on the 
idea of matching the interacting force between the streams has also been 
developed for quick estimation of Mach stem height in overexpanded nozzle 
flows. This simple analysis seems to yield reasonable results even out- 
side of its applicable flow regime. 
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NOMENCLATURE 


a.,b ,c ,d 

A jB ,C jD, 
E,F»G 




«2‘«3a’«4a 

I 
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M 
M* 
MS 
h - 


N 

P 

q 

R 


coefficients defined in Eq. (39) 

functions j see APPENDIX A and Eqs. (24) and (25) 

functions, see APPENDIX B and Eqs. (35) and (36) 
specific heat at constant volume 
distances defined in Fig. 33 
incident shock 

reference length (here maximun body height above the 
centerline of symmetry) 

Mach number 

Mach number based on speed of sound at critical condition 
Mach shock 

dimensionless coordinate normal to the body surface 
measured from the surface 

number of divisions on the starting characteristic line 
dimensionless pressure 
dimensionless magnitude of velocity 

reflected shock, dimensionless radius of curvature of body 
surface 


r 


dimensionless radial coordinaxe in the axisymmetric method 
of characteristics 


s 


SL 


dimensionless coordinate along the body surface measured 
from the shock attached point ; entropy per unit mass 

slipline 


T 


triple point 


u,v 


V 

n 


,v 


dimensionless velocity along horizontal and vertical 
direction in internal flow, see Fig. 13b 

dimensionless velocity normal and parallel to. the body 
surface 



V 


V. ■ 
V 

X,Y 


x,y 


Y 

e 

Y 

o 

z 

Y 
K 

(~) 


velocity vector, 
dimensionless velocity 

dimensionless horizontal and vertical coordinate for 
flow inside a nozzle, see Fig. 13b 

dimensionless horizontal and vertical coordinate for 
supersonic flow over an arbitrary body, see Figs. 12 
and 13a 

dimensionless nozzle exit height above centerline of 
symmetry 

dimensionless Mach shock height above centerline of 
symmetry 

parameter defined in Eq, (32) 
del vector operator 

dimensionless shock radius of curvature 
dimensional quantity 


GREEK LETTERS 


a 

Y 

6 

e 


Mach angle 

ratio of specific heats (y = l.*^ here for air) 
deflection angle 

dimensionless shock layer thickness measured normal to 
the body surface 


e 


angle from axis of symmerry to surface tangent; flow 
angle in the method of characterisrics 


K 

X 

p 

O 

CO 


angle from normal to centerline to shock tangent 

angle defined in Eq. (51) 

dimensionless density 

shock wave angle 

P/(p"^), see APPENDIX A 

Prandtl- Meyer function; Crocco number = q/c ,, 
APPENDIX A 


see 
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SUBSCRIPTS 
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1,2,3,4 


a 
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e 
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max 

Md 

MS 


n 


n3 ,114 


P-H 

r 

rS, 

s 

s 

o 

sk 

sZ 

t 

u 

w 

I, II 
e 


stagnation sxate; starting place at triple point; 
condition on body surface in the method of integral 
relations 

region 1, 2, 3, and 4, see Fig. 12; integration step 

freestream condition 

ambient condition 

centerline of symmetry 

nozzle exit 

free- jet boundary 

incident wave; integration step index 

lower part of si ip line 

maximum state 

Mach' disc 

Mach shock 

normal to streamline- coordinate 

extreme cases of cr^ = 90 degrees or 0^ = 90 degrees at 
triple point 

Prandtl- Meyer expansion 

reflected w’ave; reference state 

limit of regular reflection 

streamline direction 

body surface 

imbedded shock v;ave 

sonic condition after shock 

throat 

upper side of slipline 
wall condition 

characteristics of family I or famiily II 
conditions just behind the shock wave 
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1. INTRODUCTION 

The study of supersonic flow necessarily involves the consider- 
ations of normal or oblique shocR waves since almost any irregulari- 
ties will produce one or more shock waves in high speed flow field. 
Although the appearance of these waves can be e>:plained and taken into 
account in many of the flow situations , one of the unsatisfactory and 
yet unanswered areas involving their reflection and mutual interaction 
with the flow field is "Mach Reflection" or irregular reflection. This 
type of phenomenon was originally observed by E. Mach [l]* and later 
brought to public attention by von Neumann [2]. The role of shock wave 
interactions in the derermination of the overall flow field can perhaps 
be best appreciated by quoting von Neumann’s remark: 

"No 'situation Involving t^-ro o'r more shocks can be fully understood 
without solving such problem., and a majority of .the practically 
or theoretically significant situations involve several shocks." 
When an oblique shock I is generated within a supersonic flow 
field, such as shovm'in Fig', la, it must be reflected from the lower 
wa.ii, and the flow behind the incident shock must undergo a change of 
streamline angle 6 after this reflection. If this angle 6 is less than 
6 corresponding to the flow Mach number M_, the reflection of 'the in- 

iiiaX . • /. 

cident wave is regular — both the incident and the reflected v?aves I • and 
R are straight and intersect with the solid wall. On the other 'hand, if 
6 is larger than the maximum allowable turning angle for the flow Mach 
number a completely different situation appears (see Fig. lb). Some- 
where along the incident wave a triple point T occurs . A reflected 

^Numbers in brackets refer to entries in REFERENCES. 
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curved shock E. and a strong curved Mach shock MS (Mach stem in ti'jo- 
dimensional problem) appear at this triple point. The flow behind the 
Mach shock is subsonic; the flow behind the shock vrave R may be sub- 
sonic or supersonic. A slipline also emanates from this point indicat- 
ing"the discontinuity in entropy. This type of shock pattern is named 
Mach Reflection; its configuration also fully illustrates the compli- 
cation of the flow field arising from the multiple interaction of 
curved shocks . It should be mentioned that the incident shock remains , 
nevertheless, to be straight for the case shown in Fig. lb for uniform 
incident flow, 

Mach reflection also occurs 'within the overexpanded and under- 
ejqjanded 'two-dimensional nozzle flow field. As it is well knovm that 
for slightly overexpanded nozzle flov?, an oblique shock which is gener- 
ated from, the comer of the nozzle will be regularly reflected from the 
centerline (Fig. 2a). If the nozzle is highly overexpanded, i.e. higher 
back pressure, the shock can no longer be regularly reflected and a Mach 
reflection pattern such as shown in Fig. 2b appears. 

In both overexpanded and underexpanded axisymmetric nozzle flovis 
the irregular reflection patterns occur as well and have appearances 
essentially similar to those of the 'two-dimensional counterparxs . Re- 
ferring to Fig. 3 (an underexpanded axisymmetric nozzle): If the back 

pressure is much lower than the jet pressure at the nozzle exit, the 
compression -waves resulting from the reflection of expansion vxaves from 
free jet boundary coalesce into an incident shock v?ave inside the jet 
and the Mach disc eventually appears. 

Mach reflection also takes place in many other flow situations 
such as the transient flow field vjxthin a shock tube, sxeady flow field 
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within diffusers and nozzles and the transient shock interaction result- 
ing from simultaneous explosion of -tvro charges. Close to a wall in 
boundary layer-shock wave interactions , Mach reflection has often been 
observed. 

The Mach reflection study finds its application in many military 
systems. The flash appearing in the issuing gas jet at the muzzle of 
a gun or at the exit of a rocket motor is often attributed to the Mach 
reflection phenomenon. The investigation of recoilless and noiseless 
rifles is closely related to the jet structure and shock interaction 
problem in supersonic flow. Since 'Mach reflection would generally re- 
s,ult in a higher pressure on the reflecting surface than regular re- 
flection » the maximization of the damage effect in conjunction with 
Mach reflection was studied and utilised in the atomic bomb explosion 


at 2'Iagasaki and at Hiroshima, Japan, during World War II [3j. Rocket 
motors, missiles, jet engines, and other propulsive devices flying "at 
very high altitudes , for example , always have Mach configuration associ- 
ated with them. ' The relaxation and radiation behind the Mach shock are 
.sp. 'intense that they cou'ld be detected by a modern electronic device 
,(e,.g. ,, an -infrared detector [4]). The thoro'ugh study of the Mach rsflec- 
-tion problem t?ould, therefore, he helpful in the search and detection of 
these flying objects. 


Problems in the design of rocket motors, jet engines or turbine 
blades cannot be fully understood without a kno’vfledge of exhaust plume 
flow field of a gas jet vjhich, in turn, relies on the study of Mach re- 
flection in high speed flow. 


Academically, the solution of Mach reflection problem is of high 
interest. It not only provides useful information concerning multiple 
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shock wave pattern, but it fills the gap existing 5n the theorexical noz- 
zle analysis when the back pressure outside the nozzle is at off-design 
condition . 

The present investigation is restricted to the inviscid study of 
Mach reflection and the associated flow field of steady t^'To-dimensional 
overexpanded and axisymuietric undersxpanded nozzle flows. The experience 
gained in the stud 3 ^ of ejector flow problems at the University of Illinois 
at Urbana-Champaign [5] suggested that the problem at hand belongs to the 
category of inviscid interaction betvreen multi-streams as long as the 
central- core flow is distinct while the viscous -effects (such as prevail- 
ing along the slipline or along the jet boundary) can at most contribute 
to a modifying influence to the flov? patxern within this flow regime. It 
is thus believed that the shock confi^ration as well as the accompanying 
flow field (such as Mach stem height or Mach disc radius) can be estab- 
lished under the consideration that the central core flov? behind the Mach 
shock should eventually assume a state which is equivalent to "choking" 
in the uniform one -dimensional flow. Also, .previous studie.s of shock 
wave-viscous layer interaction [6] established the conviction that the 
shock structure can only be modified ^;ithln the rarefied flow regime. 
Therefore , the Rankine-Hugoniot shock relation for air flow can be ap- 
plied ifith confidence for the triple point in the study of Mach reflec- 
tion. 

It is pertinent to point out that , in most of the overexpanded flow 
situations under the present investigation, the back pressure is so high 
that the flow would separate from the nozzle wall. It vjould produce a 
Mach reflection pattern v?ith non-uniform approaching flov? conditions. 

It is believed xhat these modifications of the flow field, due to the 
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viscous interaction, can be examined only after the flow models for Kach 
reflection with uniform approaching flows and their methods of analysis 
are successfully developed and established. The present investigation 
is, therefore, restricted to the inviscid flox^ field only; even the pos- 
sible influences of the viscous effects are well recognized. 
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2. LITERATUHE REVIEW 

Although Ernst Mach (1838-1916) did not bring forth theoretical 
solutions , bis observation and discovery of experimental facts in com- 
pressible fluid flow and ballistic research have been the inspiration 
and foundation for many theoretical studies in- -gas dynamics.. Among 
his various findings in physical science > the irregular reflection 
pattern of shock wave interaction was first observed by him in 1875 
and was later named "Mach reflection" in honor of his excellent experi- 
mental contribution Cl], By igniting ti-ro electric spark gaps simultane- 
ously to produce two intersecting blast waves in the physical laboratory 
at the University of Prague, Mach was able to visualize a trace of triple 
point of forked shockwaves (die gegabelte Verdichtungsstosze) on a flat 
plate covered with a thin sheet of -soot. The intersecting vrave phenome- 
non of the irregular reflection configuration in a supersonic air jet 
exhausting from a high pressure reservoir to low pressure surroundings 
was also discovered by Mach using the then already-kno'im schlieren tech- 
nique, and the first schliei'en photograph of Mach configuration from nozzle 
flow was published in 1889. 

It was much later (during World -War -II ) when von Neumann brought 
this subject into general attention. The importance of study of produc- 
tion, equilibrium, and interaction of shock ^.'■aves in the science of 
fluid dynamics v;as emphasized by von Neumann. In his expository paper ' 

[21 he elucidated the basic diffei'ence betv;een regular reflection and 
Mach reflection and questioned the effect of viscosity and heat conduc- 
tion as an essential part in the study of Mach configuration because of 
good agreements 'between the theoretical regular reflection analysis using 
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the conventional inviscid nonheat-conducting shock theory of Rankine- 
Kugoniot and the experimental results . Because of considerable dis- 
crepancy bet^'Teen theoretical analysis and experimental evidence in liach 
reflection for weak incident shock , he surmised that there might be 
some other unkno--7n. discontinuities existing in nature other than the 
well known discontinuities such as shock wave, density discontinuity 
(slipline), and Prandtl-Meyer - angular expansion and compression in the 
ideal fluid analysis. 

The early studies of regular -reflection and Mach reflection of shoc?-cs 
were carried out by Bleakney and Taub. They \rere mainly concerned v.’ith 
transient flovr situations within the shock tube. Excellent agreement v:as 
obtained between the theoretical and experimental 'results for regular re- 
flecxion of shocks [7]. On the other hand, Mach reflection has been ob- 
served for weak incident shock strengths and wave angles where the xhrse- 
shock theory for the triple point predicts no non-trivial solurion. How- 
ever, experimental results did support the three-shock theory’’ when the 
incident wave of Kach reflection is strong. 

The Mach configuration is also viewed by Taub [8] as shocks with 
singularity at triple point v.’here the tangent to shock is discontinuous. 

He proved rigorously that the locus behind the singular point on shock 
should be a slipline in inviscid analysis . He also, gave the expressions 
for curvature of the shock and the derivative of pressure behind the 
shock at triple point singularity in terms of flow variables along the 
shock so that the geometric properties of the shock configuration near 
singular point for the- pseudo- stationary flow could be determined. 

Fletcher, et al. , [9] revievmd the case of pseudo-stationaicy Mach 
reflection of a nearly glancing incident shock using linearized theory. 
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With their own interferometrical experimental results they showed the 
adequacy of the linearized theory in handling the pseudo-stationai^'- 
Mach reflection of shock waves at nearly glancing incidence. 

• , By using the method of power series expansion of the flow variables, 
Clutterhan and Taub [10] gave the numerical results of flow variables 
near triple point in table form and used them to determine the trace of 
the triple point (called X'Torld line) and to construct approximately the 
Madh configuration in the neighborhood of the triple point for pseudo- 
stationary flow. 

In response to the above mentioned discrepancy betr^een inviscid 
triple shock theory and experiment for weak incident shock, Sternberg 
[11] took the effect of viscosity into consideration and introduced a 
non-Rankine-Hugoniot shock wave zone (incomplete shock wave) at the 
intersection of the three shock waves based on the reason that for vreak 
incident shock the wave was thicker and that this finite shock zone was 
able to support some difference in pressure and flov; angle. The con- 
ditions of equality of pressure and of flow angle across the slipline 
were, then abandoned. With the aid of experimental results from a shock 
tube, he found that the height of this tiTO-dimensional non-Rankine- 
Hugoniot shock zone for weak incident shock v;as several times greater 
than the thickness of Rankine-Kugoniot shock wave. Using the results 
foiaid from Busemann electric tank analogy for two-dimensional subsonic 
flow in distorted hodograph plane, he discovered the reflected and Mach 
shocks v?ere so strongly curved at the triple point that the shock wave 
angles beti-;een theory and experiment could not be compared accurately 
because the region involved was too small to be visible by the present 
experimental method. The introduction of non-Rankine-Hugoniot shock 
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wave zone necessarily seems to make theoretical analysis insurmountable 
for weak Incident shock vrave. 

The supersonic gas jet has been studied by many researchers after 
World War II. The formation of shock ;:ave in a two-dimensional super- 
sonic gas jet was examined by Pack [12] using the method of characteris- 
tics. The origin and position of shoc’c wave imbedded inside a free jet 
boundary for underexpanded nozzle flow vrere located based on the reason 
that the trace of shock was the result of intersection of compression 
waves of the same family. He observed that for tv/o- dimensional gas jet, 
the regular reflection still existed at rather high pressure ratio where 
Mach reflection pattern already appeared in the axisymraetric jet. How- 
ever, he did not go further to treat the irregular reflection problem. 

Some experimental investigations have been performed by Ladenburg, et al. , 
C13], using interferometric study, for the axis 3 anmetric gas jet from an 
orifice instead of a well designed nozzle . 

By introducing a one-dimensional fictitious nozzle extension, 

Adamson and Nicholls [14] employed the pressure distribution on the center- 
line of the flow behind an axlsymmetric orifice based on 'the method of 
characteristics to find the location of the first Mach shock. This Mach 
shock was assumed to exist at the end of this fictitious nozzle extension 
where the airi>ient pressure would be produced behind the normal shock. 

The analysis is very simple but it is not very helpful in the determin- 
ation and understanding of the overall flow field. Eastman and Radtke 
[40], hovrever, employed a different speculation that the location of 
the normal Mach shock wave coincided with the point of "minimum" pres- 
sure behind the imbedded shock vrave. 

The criterion of sonic condition at the throat in the interacting 
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flovf field was used earlier by Chow and Addy [5] in the study of the 
mutual Interactions between primary and secondary streams of a super- 
sonic ejector system. A similar idea has been adopred by Ashratov [15j 
to find the Mach disc radius of a jet leaving an overexpanded axisymmetric 
nozzle flow,. By using a one-dimensional approach for the entire core flow 
and applying the method of characteristics for the supersonic flow above 
the slipline, he was able to derermine, iteratively, the Kach disc radius 
which agreed well with experimental data. The theoretical results were 
given for M =-2.8 and M = 3.2. 

Love and Grigsby [16] made an extensive experimental investigation 
on the axisymmetric jets exhausting from sonic and supersonic nozzles 
into still air and into supersonic -stream. Some theoretical consider- 
ations were also given in their report but were not complete as far as 
Mach configuration study v;as concerned. Their experimental data are to 
be compared here with the numerical solutions. 

The method of integral relations plays an essential role in the ap- 
proximate solution of fundamental hydrodynamics equations . Traugott [17] 
carried out a one-strip analysis of integral method of Belotserkovskii 
[18] for supersonic flow over arbitrary blunt axisymmetric bodies. The 
higher order approximation using the gradient method to evaluate the flow 
variables across the shock layer was also .mentioned and'v/as different from 
the higher approximation scheme used by Belotserkovskii. His one-strip 
analysis [17] is to be modified and applied in This study for solving 
supersonic flow over arbitrary teo-dimensional body with attached shock 

k 

wave . 

South and Revjman [19] used the merhod of integral relations to solve 
real gas flow problems . They found that this method was generally ap- 
plicable to supersonic flow past a pointed body \fith vreak attached shock 
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wave. The effect of the surface flow properties propagating noriT>al to 
the boundary surface rather than propagating along characneristic curve 
was found to be self-compensated in the method of integral relations. 

They also gave the stability' criterion for determining the step size of 
integration in hyperbolic domain. 

The method of integral relations has also been adopted to solve 
the nozzle flow problem. Liddle and Archer [20] gave the equations for 
one- and two-strip analysis derived from continuity, energy, and irro- 
tationality relations. Hewlett [21] used the one-sxrip analj>-sis originat- 
_ing :^om continuity , energy , and normal momentum equations to investigate 
nozzle flovrs. His analysis is utilized for solving tV7o-dimensional 
central core flow for the present study. 



12 


3. 'TRIPLE POINT AJID FLOW RBGINE FOR NACH REFLECTION 


The Mach configuration is characterized by the rautual interaction 
of several shocV waves . This distinguished feature is manifested 
through the appearance of a triple point where shocks of different 
strengths intersect each other. The flow^ condition in the neighbor- 
hood. of this triple point has the dominant effect on the downstream ' 
flow field because of the necessary existence of subsonic flov? behind 
the Mach shock. The study of the equilibrium of flow variables near 
this triple point is, therefore, indispensable in the- Investigation 
of Mach reflection. 

For nozzles operating under the condition v;here Mach configi:iration 
occurs, a triple point T accordingly appears (see Figs, lb, 2b, 3, and ' 
4) . The shock MS belongs to strong shock family and becomes a normal 
shock on the centerline while the reflected shock R can be vjeak or 
strong and is also curved. The strengths of .reflected shock. and of 
Mach shock in the immediate vicinity of -triple point T are determined 
locally by the fact that the pressure is continuous and the flow angles 
are the same across the slipline which results from the entropy differ- 
ences in state 3 and state 4 for this inviscid analysis . Upon combin- 
ing oblique shock relations [223 > the flow regions 3 and 4 ax point T 
for an arbitrary incident shock (arbixrary y, M^, and 6^) can be solved 
from the following equations: 


2Y ,.2 . 2 ^ T - 1 

y-n S ^3 - frr 


( 1 ) 


^4 „ 2Y . 2 ^ Y-i 

P ■ Y 1 ‘l S ~ Y + 1 


( 2 ) 
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tan 6„ = 


2 2 ' 

iq sin CTg - 1 


3 tan cf„ ^ ^ 

3 M^Cy + cos 2cfg) + 2 


tan 6,. = 


2 2 

iq sin - J 


4 tan a„ ,,2, ^ . 

4 Hj^Cy + cos 2cr^) 4 -2 




P = P 
3 4 


(3) 


(^) 

(5) 

( 6 ) 


The solutions are found from Newton-Raphson ' s method for simul- 
taneous nonlinear algebraic equations by numerically evaluating partial 
derivatives [29 3. These equations have trivial solutions corresponding 
to M = M = M , and three shocks merge into one shock which is of no 

^ V 4" 

interest and is discarded. The nontrivial soluricns are plotted in 
Figs . 5 , 6 3 7 3 8 3 and 9 . 

It is interesting to note that the state behind the reflected 
shock may be subsonic or supersonic near the triple point. The con- 
dition when sonic flow behind the reflected shock occurs can be de- 
termined, in addition to Eqs. (1) through (6), for any given and y 
from the following equations; • 


tan 0^ = v 

2. tan ,,2 


(7) 


2 Mj^Cy + cos 20^) + 2 


(Y + 1)^ sin^ ^2 " ( ‘'^l ^2 


,2 . 2 


= 


|^2YmJ sin^02 - (y - D (y ~ D f'' 


2.2 

sin 0^ + 2 


> 


( 8 ) 


.2 if 

sin 0 = 2 i ' 

4yM, I 


(Y + 1) ^2 - (3 - y) 


+y<Y + 1) [Cy t 1) - 2(3 - y) f (Y + 9) I 


(9) 



These soluxions sr>e also .given in ,Figs . 5 , 5 , 7 , 8 , and 9 and ane 
labeled H =1. 

T;^o extreme cases are of interest: One is sin = 1.0, which im- 
plies that Kach shock is a norr:al shock' at triple point; that is, 

^2 = ^ 2 ’ ^4 ~ ^ pseudo-regular- type reflection as shovrn in 

Fig. 4a. The other case is sin c?^ - 1.0, which indicates that the 
shock between state 2 and sxate 3 is a normal shock (see Fig. 4b). 

This is the physical limit U'here a solution to the triple poinx exists 


under the condition of fixed and Y* 


The solutions to these tv?o extreme circumstances are easily 
found from solving Eqs.. (l) through (8) with additional restriction of 
sin = 1.0, 6^ = 0, or sin c,g = 1.0, 6^ = 0 for the specified value 
of Y and M., and are shown in Figs. 5, 6, 7, 8, and 9 identified i-rith 
= 90 degi'^ees and cr^ “ degrees, respectively. It should be noted 
that for < 1.434 ^^(y t 3,)/2 as given in [42]^ with y = 1.4, 'no 
triple point solution is possible for this nozzle flow, and for 
M.^< 2.406, M.., will be subsonic regardless of the value of 6 . 

The numerical values of triple-point solutions together with 
their two extreme cases for various i^ach numbers are tabulated in 


Table 1 through Table 13 v;here 1, 2, 3, and 4 refer to the different 
states at triple point as shown in Fig. 4, DELTA is the deflection 
angle in degrees , SIGMA is the wave angle in degrees , M is the Hach 
number, P2/F01 is the ratio of pressure at state 2 to the stagnation 
pressure at state 1, and (P4/P1)* is the ratio of the pressure at 
state 4, v.'hen expanding to sonic condition, to the pressure at state 1. 

In order to study the f loi? field associated vfith Mach reflection , 
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it is 'also necessarj.' to identify the i''egirne where Mach reflection occurs 
It is well known, according to [7], that for an incident shock wave 
striking a horizontal lower vrall (see Fig. la) , regular reflection vjould 
occur if the turning angle- 6 associated with the incident shock is less 
than sonie limiting value The relations- bet^/een flow Mach nuirber *■’, 

6^^, -and wave angles are 


tan 6 


2 2 

M, sin a . - 1 
1 1 


vl tan O. , ,2 f \ 

1 M.J Cy + cos 20^) -r 2 


tan '6 


. 2 

„ ,sin a - 1 

2 • -2 r 


rit 


uan cr^ M^(y -r cos 2o ) + 2 
2 r • 


( 10 ) 


( 11 ) 




(Y + 1)^ ~ *^i **■ 

2YM^ sin^ a. - (y - 1)1 F(y - D sin” cJ.'t 2 

X X XX 


sin O = 
r 


. 1 


4ym; 


r 

{ 


(Y + 1) iC - 4 


( 12 ) 


(13) 

The solutions are given in Fig. 10a for y = 

Since must be greater than one in order that there may be a 
reflected shock , the upper limit is thus tentatively established under 

the condition of the sonic state behind the incident shock and is 

♦ 

marked as curve 6 . in Fig. 10a. 
s 

If the two extreme results, 6 „ and d ,,,.from triple-Doint so- 

’ n3 n4’ 

lution are also plotted in Fig. 10a , there is some overlapping area 

belov; curve 6 „ vrhere both regular and Mach reflections are possible 

and some empty area between curve 5 . and curve 6 „ v;here regular 
^ sx. n3 


+ 1)|^(Y t 1) + 3 (y - 1) + 16 j- 
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reflection is iisnossible . The enrory area beri^een curve 6 „ and curve o 
above regular reflection limit, curve is > perhaps, due to the limit- 

ation of the inviscid assumption. For this lo-^er Mach nuniber regime 

- 1.0), Sternberg [11] discarded xhe criterion of equality of pres- 
sure and of flo\? angle acr-oss the slipline and introduced a non- 
Rankine-Hugoniot shock at the triple .point Trirh “he viscous effecx con- 
sidered. From an electric tank experimenr , he ’.^as able .to find the so- 
lution to the Mach reflection problem. Moreover, holder [231 proposed 
a smooth curving shock pattern for the case uhere neither regular nor 
Kach reflection Is theoretically possible. 

In the overlapping portion both tvro-shock and three-shock theories 


predict nontrivial solutions. Whether fwo-shock theory or three-shock 
theory should be used depends, at ‘least, on the boundary conditions. 

As indicated in Fi.g. 11, if a plug vrith a particular shape shotm vjere 
. inserted at the exit of a two-dimensional nozzle giving a two-shock 
theory pattern, the removal of this very plug might arouse a Mach con- 
figuration from dormancy under the same ambient pressure. Hypothetically 
this plug could be used to control the appearance of Mach configuration. 
For the problem shown ;.n Fig. lb, the experimental evidence of Holder 
[23] shows that at M., = 2.8, the transition from regular to Mach reflec- 
tion occurs at the limiting case of three-shock theory rather than xhat 
of tiio-shock theory and that the Mach shock is a normal shock during 
this transition. Therefore, the stead 3 '‘ Mach reflection produced from 
nozzle flo'w of this study is confined to the domain bet^.'een curve 6 

and curve 6 . in the nresent analysis. 

n4 ■ •’ 

The relation between je~ deflection angle at nozzle exit 6 and 
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the ratio of anibieiit pressure-to-nozzle stagnation pressui’e 
given in Fig. 10b for various- Mach 'nuifbers. For convenience of refer- 
ence, the region specified in Fig. 10a is also plotted here. At high 
Mach nurrbers the flo;-? is relatively? sensitive in the sense that a very 
small change in ambient pressure results in tremendous variation in 
the jet deflectibn angle 6^ at the nozzle exit. 
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4 . FUNDAMHNTAL EQUAT10j:S 

For* a steady inviscid gas flow, xzithout shock discontinuity inside 
the T>egion considered, the conservation equations in vector fom are: 

Continuity : 

7 ‘(pV) = 0 , (14) 

MoiaentuTn : 

p(V • V) V = -VP (15) 

Energy : 

(V • V)f^1 = 0 (16) 

Lp^'J 

where the equation of state for ideal gas has been tacitly used and (~) 
indicates dimensional quantity. 

The first two equations simply mean that mass and momentum are con- 
served, The last equation states the entropy is constant along the 
streamline under the above assumptions . Of course the entropy can change 
across the streamlines in the case of flow behind' curved shock wavs . 

The vector form of the fundamental equations shown above' give concise 
expression for the conservation principle vfithout referring to an;/ particu- 
lar coordinate system. For different physical problepis , the conservation 
equations have different forms for specific coordinate systems used ’which 
can be derived from the vector form through vector analysis and coordinate 
transformation. 

The fluid considered in this study is assumed both thermally and 
calorically perfect with y - 1-4. 




5.i BASIC DESCRIPTION OF THE FLON' FIELD 

As sho^NTi in Fig. 12 for the case of overexpanded two-dimensional 
nozzle flow with a triple point shock sysxem iirhedced in the flo:; field, 
regions 3 and 4 are separated by r’ne slip line of a yet unknora configu- 
ration.' It is obvious tliat the ti'jo streams have to coexist so that 
along their common boundary they shall have the same pressure and flow 
direction. Reglbn 3 is usually supersonic (although it may be subsonic 
close to the triple point T) w'hile region 4 is invariably subsonic. It 
resembles a great deal the problem occurring within a two-dimensional 
ejector system [5] . It appears to be extremely attractive to approxi- 
mate the lov?er stream b 3 ’- a one-dimensional treatment » especially in view 
of the fact that the variation of entropy seems to be insignificant. 
However >5 it has been observed that in the early part of the flcu' develop - 
men'l: the pressure is increasing along the slip line. A one- dimensional 
flow analysis for the lower stream will not be able to match such a flow 
condition, and a more sophisticated flow analysis is thus needed in this 
region near triple point. 

The flow field of region 3 for the upper stream is bounded by 
a curved, attached shock and is to be treated as a supersonic flow past 
a curved body whose configuration is still unknown (see Fig. 13a). It 
presents itself as a part of the flovj field produced by a supersonic 
flovr past a pointed body, .Although the method of characteristics will 
be useful when the flow is supersonic, it v;as originally expected rha'c 
the genera] method of anaij-sis for blunt body problems will be needed if 
the flovr 'is subsonic behind the curved shock. The mexhod of integral 
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relations developed fay BorodnitsTO and applied fay Eelotserkovsfaii [121 . 
to fluid flovj- profalenis is thus adopted for faoth the upper and lo’.ver flo'..’ 
fields -during the early part of the flow development. 

From the fundamenral equations uhe momentum equation along the 
streamline direction can be integrated to obtain the Bernoulli's equation. 
Due to different coordinate systems used, one has different sets of 
equations for external and for internal flow. Hence, one obtains the 
follovring two-dimensional equations in streamline coordinate system: 
External Flow for Region 3 (Fig. 13a) 

Continuity : 


3n 



R t n 
R 



0 


(17) 


Normal Momentum: 


9 

R -r n 1 

3 , — ~ . 

ov 
• s 

3 I R -r n 



•r — -■ ypv V ; 

— — — — “ 


3n 

L " R J 

8S ® 

R 

3n L R - 


NOTE: These'-are Tollraien's equations [243 for 

flow. 

Internal Flow for Region 4 (Fig. 13b) 

Continuity'- : 

^ ^ = 0 
3x Dy 

Normal Momentum. : 

— ov — ■ 8v 3? 

pu — + pv — + ■ — = 0 

3x 3y 3y 


(18) 

inviscid 


(19) 


( 20 ) 


The Bernoulli relationship can be conveniently expressed by 

-5- ~ q^ = constant along the streamline (21) 



and the entropy equation is 


D_ 

Dt 



= 0, where’ 



the substantial derivative. 


( 22 ) 


-Sqs. (21) and (22) are used for both external and internal flows. Here, 
C^) refers to dimensional quantifies. 


5.2 THE METHOD OF INTEGRAL R5L.ATI0NS 

The basic procedures of applying the method of integral relations 
[17,18,21,25] fo this study are to express the governing nonlinear con- 
tinuity and momentum equations in streamline coordinates , to \jrite them 
in divergence forms , and then to integrate them along the normal to 
streamline direction. The unknown integrands are approximated by some 
polynomial functions (or some other convenient continuous functions). 

With the aid of Leibnitz's rule for differentiation under the integral 
sign, the original nonlinear partial dlfferenrial equations are re- 
duced to a set of ordinary’' differential equations readily solvable by us- 
ing a digital computer. Since the integration is a smoothing process, . 
the overall character of the solution is preserved even though the de- 
tailed description of the flo^.' field cannot be achieved. The accuracy 
of the solution depends on the number of strips of the normal direcrion 
divided . 


5.3 SUPERSOITIC FLOW OVER TWO-DIl-ENSIONAL ARBITRARY BODY 

If the normalizing -DaraTneters a , ? , p , ar:d L are intro- 

'IT3X Ooo Oco 

duced to normalize the flow parameters in a one-strip analysis , and if 
the unknora integrands are approximated by linear function in normal 
direction 

f(n,s) = f^(n,s) + ~ [f^(n,s) - f^(n,s)] 


(23) 



the nonlineax* partial differential equations, Eqs, (17) and (IB), for 
external flow can be integrated along the normal direction from the 
body surface to shock layer thickness e and, after tedious mathemati- 


cal manipulation, will yield: 


G + E 

d.< _ cs 

ds F 


from normal momentum equation. 



from continuity equation , and 



(24) 


(25) 


d£ 

ds 


1 + 


R 


tan (k + 6) 


(26) 


which is a shock geometric relationship, see Fig. 13a, 

• Functions A, B, C, D, E, F, and G are complicated functions of 

K, , E, Y> 0> d6/ds , and are given in APPENDIX A. Three un- 
o 

known s , fC, V , and e, are left to be determined. It is interesting 
o 

to note that the above t>?o- dimensional relations can easily be deduced 
from Traugott's axisymrr.etric equations [17] by allowing the body ra- 
dius to approach infinity. 


5,3.1 Detached Shock Nave 

If the shock wave is detached , the initial 
the above coupled first-crier nonlinear ordinary diffe: 
are ' . 


ccnai 
'ent ' a 


on 

ec 


C"’~ 


K 


Is = 0 = ''s 


= 0, . e 


s = 0 “o 


(27; 


s = 0 


The shock stand-off distance e is unknov:n and must be solved as cart 

o 
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of the solution. There are tx-fo singular points involved. One occurs 

at B = 0 and the other at F = 0. The singular poinr of B = 0 occurs 

ar = Jiy ~1 )/(y t~l)' v;hich is the sonic condition. Ihis singular 

behavior has been adequately used as the criterion for iteratively 

determining the correct stand-off distance £ for axisymmetric bodv 

o 

[17,18,25], since the surface velocity gradient must be finite for 
the smooth body profile. Cn the other, hand, sone investigators [25] 
have questioned the applicability of the one-strip method of integral 
relations to tv'o-dimensional detached shock wave problems due .to The 
ill behavior of the second singular point F = 0. This wdll be dis- 
cussed later. 


5,3.2 Attached Shock Wave 

In this study of Mach reflection , the shock vrave is al- 
ways attached to the triple point,, A .limiting case corresponding to 
the zero shock stand-off disuance must be investigated, since the fore- 


going equations assume undetemsined form if xhe condition e = 0 is im- 
posed direcxly .into these equations. By formally using L’ Hospital's 
rule and by noting v 0 , P P at the attached point when £ 0 , 

the coupled differential ecuaricns , F,‘s. (2*0, (25), and (2G), becc v. 


[£], 


tan (0 T k) 


2 2 

V tan (0 + k) - 
s 

o 


Y -i - 1 2 

1 - V 

. s 

Y - 1 o 
1 - v^ 


V k k„ 
S 1 2 


tan (6+K)r(l-^v^)(^+-^^-4 )+v^ k, k^ tan (Stk] 

. . - o o 1 - V o _l 
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with 


* 1 ■ 

R 


■ dK * 


d8 




- ds _ 

e -»■ 0 

- ds - 


sin (k + 0),froin shock geometry 


(30) 


e 0 


rdv 


ds 


0) 


TO ds 


r, \ ■''s 

“ t— 2— 


2Y<j._ 


1 - V 




tan (6 T k) + V k. k- 
s o 

o 


^^-0 u 


V k^ k^ 
1 3 

o 


(31) 




The parameters have the same meaning as those of APPENDIX A. 

These are the equations for the initial derivatives of the unknowns 

e, K, and v The initial values of K and v^ can he found from obliaue 
’ ’ s ®o 

o 

shock relations. These equations can also be iiiterpreted as the relatioii- 

d8 


sur- 


ships bet^'Teen shock curvature 1/R and body curvature ~ Hi" ) 

-face velocity gradient dVg^/ds (or surface pressure gradient) 'at the ar- 
tached poinu of shock wave on the body surface. Figures 14, 15, 16, and 
17 present the behavior of these relationships both for weak and for 
strong shock waves. A. discussion on the general strong shock relations 
can also be .found in Reference [26]. It is noteworthy that the exact 
expressions for weak shock wave have been derived by Thomas [27], Taub 
[8],. and by Lin and Rubinov [28] from different methods of analysis for 
uniform approaching flow, and by -Holder [23] for non-uniform incident 
stream with rotationalitv in front of curved shock viave. 


5.3.3 Treatment of Singular Points 


If one sets 
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For a fixed value of z > 0 , there corresponds tv?o different 

positive values of v : one for subsonic flow and the other for super- 

®o 

sonic flov?. At sonic condition, ( -'(dv ) = 0. Therefore, near sonic 

3 

' o 

speed, "for a knovm value of z, the .responding v cannot be determined 

o 

efficiently from Wewton-Raphson' s ; ethod for solution of nonlinear alge- 
braic equations. Other, methods such as the bisect method or the false 
position method [29] vrill yield the solution near sonic point. At sonic 

condition, for y ~ 1.4, z = 0.2721655 when v = 0.408248. 

®o 



I'he ill behavior of the differsatial equation and the inaccuracy 
of interpolating too close or too far near sonic condition as mentioned 
in [17,253 have been eliminated. The singular behavior of the differ- 
ential equarion has changed, as a resui t of 'the transformation shcun 

I 

above, to the extreme behavior of xhe nonlinear algebraic equation 

which is much easier to handle. Also, because a small change in s -will 

cause tremendous change in v near sonic poinx, a' smaller integration 

step size should be used there. Since z is an odd function of v , 

there is another branch of solutions giving negative values of v ^;hich 

o 

is physically unrealistic and is ignored. 

By noting that in teo- dimensional flovr F/e is a function of ic only 
V7hen Yj and 6 are independent of K, it also seems to be possible to 
use the similar treatr.enx as before to eliminate the singula-rity occur- 
ring in differeriXial equaxions ,v7hen F = 0. nowever, because of the ap- 
pearance of dK/ds in 5q. (25), nothing will be gained when the technique 
utilized above -is applied to Eq. (24) — even the complexity of the func- 
tion F/e is ignored. Nevertheless, since F/e = f(ic), if y, and 6 

were held fixed, one can investigate, for e 0 , i?hen F/e ’.fill vanish. 

Figui-’s 19 is the result of the solutions of F/e = 0 for y = 1.4. 

Two curves from cbliaus shock ecuations for M =2.0 and H =5.0 are 

* “ 07 CO 

also plotted in the same figure for comparison. The intersection of 
these curves with solution curves of F/e = 0 at the same ficv? Mach 
number occurs at the attachment limit which is a particular soiucion of 
F/e = 0 and is also one of the starting ppints for attached shock v;avt- 
problem. Since shock is a local phenomenon., the magnitude of the ve- 
locity vector ju$t behind the shock wave at point awaj? from body is 
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independent of the surface tangent angle 0, that is, = f(H^,K) only. 
Therefore, the flow properties just behind the shock wave can be found 
by tracing oblique shock polar for the corresponding free stream hach 
number. Noting that the sonic state is very close to the attachment 
limit point on 9-iC diagram for the same free stream Mach number, one 
has, in general, t*wo different categories, depending on xjhether the 
state behind the shock is supersonic or mixed type. 

1. Supersonic State Downstream of Shock Wave 

Because K will not decrease in general for the unbounded flov; 
field and the maximum value of tc is the one when shock weakens 
to become a Mach vrave, F/e will not vanish for different com- 
binations of 6 and tc along shock polar as illustrated in 
Fig. ISa. 

2. Subsonic and Supersonic Mixed State Downstream of Shock Wave 
As shown in Fig. 19b, the singular point F/e = 0 will be un- 
avoidably met for all the smoothly varying body profiles carry- 

. ing flovf changing from subsonic to supersonic state . 

• From these results the conclusion is xhat the one-strip method of 
integral relations is not suitable for solving t-vro- dimensional blunx 
body wixb 'detached shock or pointed body with strong attached shock. 

The two-strip or multi-strip method is required and the calculation 
will be much more involved. This agrees with the comments made in [25]. 
However, for solving two-dimensional pointed body with weak attached 
shock, the one-strip method of integral relations is generally applicable 
(see also [193) for slowly changing body profile v;ith nearly constant 
value of K such as the case in our study of Mach reflection. 

In conclusion, in the one-strip method of integral relations for a 
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■two-dimensional profile the singular point B = 0 occurs at The sonic 
condition and can be handled i^ithout difficulty, but the occurrence 
of the other singular point F = 0 will greatly restrict the applic- 
ability of the one-strip analysis. 

It does not appear to be necessary to use the multi-strip analy- 
sis of the method of integral relations with the present method of 
approach in order to avoid the erratic behavior of dK/ds in the one- 
strip analysis. In weak incident shock regime, where sribsonic flow 
behind the reflected shock prevails , it is believed that the triple 
point solutions do not provide adequate initial conditions as they do 
not agree vrith the experimental results according to [7!]. Worst of 
all, the three-shock theory does not give any nontrivial solution 
while. Kach configuration has been observed in actual flo',/ situations 
with weak incident shock. If The multi-strip analysis is to be use- 
ful for weak incident shock vrave, instead of the results from the present 
inviscid three-shock theory, the initial conditions at triple point 
must be obtained from other more sophisticated approaches, for example, 
Sternberg's .non-Rankine-Hugoniot model. 

5,4 FLOW INSIDE A TWO-DII-ENSIONAL NOZZLE BY THE METHOD OF INTEGRAL 

RELATIONS 

As mentioned previously, -a one-dimensional treatment of lov;er 
part of the slipline is not adequate for solving the lo/zer fiovr field 
close to triple point. A more detailed approach is indispensable. Re- 
ferring to Fig. 13b, upon introducing the assumptions that entropy dif- 
ference along the Mach shock is negligible and after the variables are 
non-dimens ionalized by reference quantities, xhe integrands can be 
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2 ’ 2 ^ 2 
pv = (pv ) n 

W 

where n = y/y • 

The integration of the system of equations for internal flow vrill 
yield two ordinary differential equations as the consequence of the ap- 
plication of the one-strip method of integral relations. They are; 


dx - °-w - - B, 


(35) 




Aj Fj - ^1 Cl 
^1 ^1 “ ®1 ^1 


(36) 


with initial conditions . ■ 

^ ' • (35a) 

o 

q = Q' • (36a) 

^c 'C 

Q 

where A^ , , C , D, , E^ , and are functions of Y , Y' , Y", q , and 

q^, and are given in APPENDIX B. The detailed derivations and dis- 
cussions on singular behavior of these differential equations are 
given in [21]. Note that Eqs, (35) and (36) are functions of wall 
height, slope, and curvature. 

The rotationality behind this strong Mach shock is usually very 
small and can be neglected, as can be seen from the results of triple 
point solution at = 3.0, for example, for 

• 6^ = 25 degrees = 81.434 degrees = 0.551 

6^ = 30 degrees 0,^^ = 86.24 degrees = 0.491 
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in comparison to the normal shock ar centerline C = 90 degrees, 

M = 0.4752. An extension of the analysis to include the vorticity 
can easily be made and was done '[30]. 


,5.-5 THE .METHOD OF CHARACTERISTICS 


In the supersonic flow" field the method of characteristics is a' 
powerful tool. If the rotational! ty behind the curved' shock wave is 
also considered, the fundamental equations can be solved by tracing 
the physical and the hodographical characteristics simultaneously in 
a step-by-step manner. For the tv/o-dimensionai case, the characteris- 
tic equations are (see [31,32]): 

= -tan ( 8 ■+ a) 


dx 


I, II 


1 dM»'» 

de 


= + tan a - 


I, II 


. 2 ‘ d(s/c ) 
sin a V 

Y (Y - 1) d0 


(37) 

(38) 


I, II 


The last term in the second .equation accounts for the vorticity 
behind curved shock. V’hen derivatives are approximated by difference 
forms , these equations can be solved simultaneously to give the unknoxTi 
quantities and, therefore, the whole flow field (see APPENDIX C). 

Because of- the erratic behavior of dtc/ds of Eq.‘- (24), the one- 
strip method of integral relations for two-dimensional blunt body prob- 
lems is not suitable for the calculation of sxibsonic flow behind the 
triple point in region 3 , which is contrail to the original intention 
of adopting this method. Although the method of characteristics could 
be used for all of the supersonic flow computations in region 3 , the 
requirement of information' on height,, slope, and -curvature of v?all for 
quasi 'tn-ro-dimensional analysis in the internal flow renders Eqs.. (24), 
(25), and (26) of the method of integral relations to be the natural 
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choice to match with the upper flow based on Eqs. (35) and (36) at the 
early part of the calculation. The method of characteristics » however > 
will be employed for the upper flow in conjunction vrith the one- 
dimensional treatment of central core flow after the early complications 
in the core flow have been resolved. 


5.6 THE CALCULATION 'OF INTERACTING FLOW FIELD 

With the results at the triple point and rhe relations given by 
Eqs. (28), (29), and (31) as the initial conditions for exxernal flow 
and Eqs. (35a) and (36a) as the initial conditions for internal flow, 
the system of differential equations, Eqs. ‘(24), (25), (26), (35), and 
(36), can be integrated numerically. The fourth-order Runge-Kutra 
numerical method is used in the integration. For each -step of inxe- 
gration, one assumes that the slipline has the profile of a polynomial 
of third degree. 

• Y = a. + b. X + c. + d. ' ' (39) 

1X1 1 

where the coefficients a. , b^ , c. , and d. are to be determined for each 

i’ 1 i’ 1 

step i (see Fig. 12). 

For the first step of integration , a slipline curvature -l/R^ is as- 
sumed; This amounts to employing a circular arc for the first step as 
the slipline w'ith radius of curvature R^. Therefore, 

a, = Y(o) = Y 
1 o 


= Y'(o) = Y' 
1 o 
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where prima.(’) indicates the derivative taken vrith respect to x, and 


is found from solving 

20i t fl t + 2=^ Xj) 

O L 

vihen an arbitrary increment x^ is selected. One should note that 
is the Mach stem height having length of unity and is known from 
triple point solution. is yet to be determined. 

For the steps thereafter, d^ is the only unknovm since a^. , b^j^, 
and can be determined from the results of the previous step. In 
this v/ay, the slipline V70uld have no discontinuity in slope and curva- 
ture. For "an assumed R. , the unkno'.m d. can be found from 

X X 

r - 

2c. + 6d. X. =~4- 1 +fb. + 2c. X. + 3d. X?) ‘ (40) 

X 11 R. L 11/ 


a. = Y. , 
1 1-1 


b. = Y' 

1 1-1 


°i ~^2R.. 

1-1 


1 + YI , 
1-1 


Here is the x-coordinate at the end of the i th integration step and 

Yi is the Y-coordinate at the end of the i-lst integration step. 

For the selected geometric configuration of the slipline in each 

step of the calculations, integrations of Eqs. (24), (25), (26), (35), 

and (36) can be carried out. In particular, Eqs. (25) and (35) would 

give V and a and thus the pressures across the slipline at the end 
o 

of each step. The requirement of pressure matching thereby allows the 
determination of the correct R^ for that particular step (R^ for the 



33 


first step ) . 

It should be mentioned that in the first step of the calculation 
of the lower stream » the difference in the x-coordinates of the triple 
point and the .leg of the Mach shock- is ignored for the purpose of 
simplification. 

It is also obvious from Fig. 20 that the coordinate transformational 
relations between the lower and upper portions of the flow field are' 
given by . ■ ' ' 

Y = Y cos 6^ - X sin + Y (41) 

2 2 o 

X = X" cos ^2 + sin 6^ • (42) , 

The error analysis and the choice of step size in the numerical 
integration are discussed, for example, in [33^34], The simplest way 
in practice to 'obtain good results is 'to integrate the -same equations' 
twice: first using reasonably small step size, then using tv/ice larger 
step size and comparing the results. 

• -'The integz’ation process is carried out until xhe fluid along the 
slipline accelerates. Thereafter, the analysis for the lov-fer stream 
is replaced -by a ohe-dimensional isentropic flow while the method of 
characteristics is employed- for the upper flow for simplification 
pzirposes. Meanwhile, a Prandtl- Meyer expansion, which V70uld occur as 
the result of reflection of the curved shock from’ the free jet boundary, 
is inserted into the flov; calculation. In sv?itching from the method of 
integral I'elations to the method of characteristics for the upper flaw 
field, gradient method [25] can be used to recover from previous results 
of - the method of integral relations the information needed for the 
method of characteristics . In all of the calculations , the reflected 
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shock is~ observed to be rather straight. The vorticity term in Eq. (38) 
can, therefore, be neglected. 

Here one notices that another important reason for switching to 
the method of characteristics for the upper flow field is based upon 
the fact that the method of integral relations cannot handle the inter- 
section between reflected shock and free jet in a satisfying manner. 

The Prandtl-Meyer expansion zone Aw is arbitrarily divided into an 
increment of approximately one degree each where 


W(M*) = 


- /Y+A 


Y - 1 


tan 






- tan 


-1 


11 *^ - 1 


, Y - 1 x,...2 

1 — 7 ": — : — ^ 


Y - 1 ^ . Y + 1 

is the Prandtl-Meyer function resulting from direct integration of 
the two-dimensional characteristic equation, Eq. (38)., when the entropy, 
gradient term is disregarded. . ’ ■ ■ 

In the field or along the free jet boundary where pressure is con- 
stant, the unknoim flovr properties are found from conventional calcu- 
lation procedure b 3 >- tracing characteristic curves in physical and in 
hodographical planes as described in APPENDIX -C. 


The procedure .required for matching region 3 using merhod of char- 
"acteristics and region 4 using one-dimensional analysis along the slip- 
line needs special treatment. This scheme is similar to the study of 
the interaction bet^-jeen primary and secondary streams of supersonic 
ejector systems [5j. Referring to Fig. 20: If and are the di- 

mensionless reference height and 1-Sach number, and if flo;; field has 
been, calculared up to state 2 and flov; variables at states 1, 2u and 25. 
are all knovm, one may write for state 3 



■ o4- 

3 

o3 


_ L T + ^ 3u| 


[Y/(Y-1)3 


_ Y - .1 ,,,,-27 
^ Y t 1 *3i,_ 


1 / 


''Cy-D3 


(43) 


from equaliry of pressure across the slipiins, and; 

Y Q(M ) 
r r 


^3Jl 


(44) 


from consrancy of flow rate for a one-dimensional nozzle flow where 




l,{(Ytl)/C2(Y-l)]} 


(45) 


The characteristic equations, Eqs. (37) and (38), ^•fith entropy 
gradient neglected, can be written in finite difference form as 


^ - x; - = <®23> 

3u 2u 


ha -Yu n ^ 

^ <«13 - “l3> 


(46) 


(47) 


3u lu 


13 


'3u 


- 


lu 


e„ - e. 

3u J.U 


= -tan a. 


13 


(48) 


determine five unknovms X„ , Y„ , K* , 0„, and 

o\Z oVL oXl O 


One now has five equations, Eqs. (43), (44), (46), (47), and (48) to 

M - . The unknown Y, 
is related to Y^^ through coordinate transformations 'from Eqs. (41) 
and (42), and ?q4'^^o 3 ^"ound from triple-point solution. 

An early insertion of Prandtl-i’eyer expansion '.rould cause the 
lov?er flov/ passing through a minimum area with subsonic Xach number, 
while a delayed insertion vriil cause the lower flew reaching sonic 
condition at a section with decreasing area. The correct location 
vfhere the curved shock-jet boundai^ intersection occurs is such that 
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the lower stream will assuma a sonic state at a minimum area location. 
Once this condirion is met, it will be a simple matter to determine 
the location of the Mach stem height relative to the height of the noz- 
zle exit section by using triangular relations (see Fig. 12) 

sin (ISO® - - X) sin 6 




-M 




sin (O^ - ^2^ 


(49) 


Y - Y 

Y 3 ± 

e ” tan 


(50) 


here 


X = tan 


Y - Y 
P-M o 


- V 


M 


(51) 


and is the shock wave angle corresponding to and 


5.7 RESULTS DISCUSSIONS- ‘ 

The method of analysis given above is indeed fruitful for the 
present study v;ith supe 2 ?sonic state behind the reflected shock at the 
triple point. A successful calculation using double precision mode 
takes two to five minutes computation time on an IBM system 360/75 
digital computer. Figure 21 shov.’s the flow condition of the lower 
stream when the P-M expansion starts at various places for = 3.0 
and 6^ = 30 degrees . It clearly demonstrates the dependence of sonic 
state of the central core flow oh the location of Prandtl-Meyer ex- 
pansion fan. A similar trend is found for other flow conditions. 

The early occxirrence of P-M'fan, say at fig* 21, gives subson 

state, = 0.72, at the throat (o^ = 0°) for the slipline. The de- 
layed '.appearance of P-M fan, say at = 40, results in a sonic state 

at a position where the area is not a minimum, which contradicts the 
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qne-dimengional isentropic .flow theory ;'-Thle\ correct position of P-M 
fan 'is' therefore found by iteration ( approximately at Xp_j_j = 35.53 for 
the present’ case) . It is to be noticed that the iteration procedure 
is carried out only up to the state where Kach number in central core 
is close to 1.0 at the minimum area. The sonic state at throat .cannot 
be reached in an exact manner, since sonic point behaves as a saddle 
point mathematically. 

All above sfated discussions are based on the isentropic one- 
dimensional approach in central core flow after the highly non-uniform 
■state near triple point has been solved by using quasi tt/o-dimensional 
-integral method. For miore rigorous 'and accurate analysis not only the 
vorticity behind the curved Mach stem shoch should be considered in 
-the quasi t^70-dimensional integral method, but the nonisentropic effect 
in the continued calculation of the one-dimensional treatment must be 
taken into account also, In that case, the lower stream flow is not 
necessarily sonic at the throat. For one-dimensional nonisentropic 
nozzle flaw, Snyder [353, Bryant [363, and Ferguson [37] shovjed that 
-for nozzle efficiency less than one, the Mach number at throat of the 
nozzle is also less than one. On the other hand, in an irrotanional 
quasi tv7o-dimensional nozzle analysis [213, the "choking” condition 
is defined as the throat singular point which permits continuous ac- 
celeration from subsonic to supersonic flow. Ko’wever, these deviations 

are small compared with those of [21,353 for high efficiency nozzle 

♦ 

flow. 

Figures 22 and 23 give the flo;r patterns for 6^ = 30 degrees and 

= 25 degrees, respectively, at = 3.0,* For 6^ = 25 degrees the 

*There v?as a minor error in reporting the results in [38j. This has 
been corrected in Figs. 22, 23, and 24. 
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mini-Hixim area (sonic state for region in one-dimensional isentropic 
analysis) occm?s in the P-H expansion zone as sho;vh in Fig. 22. For 
^2 = 30 degress, however, the minimum area stands at the place ^?here 
the last P-H expansion wave ref leers from the free jet boundary and 
Intersects 'w'ith the slipline , since the ambient pressure is higher 
and the P-M zone is smaller (Aw = 8.84 degrees compared with Aw = 15.26 
degrees for the case of = 25 degrees) as shown in Fig. 23. 

Figure 24 illustrates the velocity distribution along the slipline. 
The decrease in velocity along the early part of tihe slipiine accounts 
for the necessity of using quasi two-dimensional analysis in region 4 
near the triple point. 

Figures 25 and 26 present the informaxion of the height and the 
position of Hach stem for tv/o-dimensional overexpanded nozzle flow. 

It is understood that at xhe same nozzle ^5ach number, it takes more 
computation to reach the sonic condition in central core flovr for lower 
ambient pressure levels. These curves could be, therefore, extrapolated 
to the zero Mach stem height corresponding to the regular reflection 
pattern. For higher ambient pressure, the Mach stem is higher and stands 
closer to the nozzle exix plane. 

The correct positions of P-M fan for various flow conditions are 
shoim in Figs. 27 and 28. It has a trend similar to that of the Mach 
stem height. Figure f- shows the range of expansion of P-M fan. 
Evidently at the same - • le Hach number lower ambient pressure induces 

wider range of P-M exp*: .*n. 

The height and Ic . ..ion of the throat of central core flow are 
plotted in Figs. 30 a ..31. The higher Mach sxem height is associaxed 
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with higher throat height for the same nozzle exhaust condition. This 
indicates that for a propuls ional device flying at a designed nozzle- 
Kach number from higher pressure atmosphere to lower pressure sur- 
roundings , the Kach stem height decreases as does the throax height 
in the central core flow. Therefore, the mass flow rate passing 
through the central core is reduced and the pumping effect of this 
central core is reduced also. However, the higher surrounding pres- 
sure does not necessarily result in a throat position closer to the 
.nozzle exit plane . As depicted in Figs . 22 and 23 , different ranges 
of P-M expansion for various amhient pressure levels would affect the 
location of the throat significantly. This perhaps accounts for the 
occurrence of kinks in the curves as shown in Fig, 31, 

Figure 32 reveals the values of M* on the upper side of slipline 
;3t the place where the lower side is at sonic' condition. It is interast- 
,ing to note that at the throat the velocity difference on the tv/o sides 
of the vortex sheet (sliplirie) decreases as the ambient pressure in- 
creases under the same nozzle design condition. 

With the data given above, it should be fairly easy to construct 
a Mach reflection pattern from an overexpanded two-dimensional nozzle 
flow, when it occurs. First of all, the nozzle exit height and Mach 
number are specified. For a particular ambient pressure one finds a 
corresponding 6^ from Fig. 10b. The properties at triple point' are 
found in Figs. 5 through 9 (or from Tables 1 through 13). Using this 

i 

information and the data shown in Figs. 25 through 32, the appropriate 
location and height of Mach stem, P-M fan, and central core throat can 


■ a] 1 be laid out . 
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5..8 A.PPR0XIMATE METHOD FOR ESTII-5ATING MACH STEM HEIGHT 

For higher back pressures , state 3 at the triple point becomes 
subsonic and the method described previously doss not yield any fruit- 
ful results, since the calculation cannot overcome the singularity 
corresponding to F = 0 (see Eq. (24))-. However, the location of the 
Mach stem and its height under these situations can be determined ap- 
proximately from a model emphasizing the balance, of. the interacting 
force bet^7een the two streams . As it turns out , this model can also 
be employed to give reasonable results even vrhen state 3 at the 
triple point is supersonic. 

One realizes tha~‘ under' this flow situation, the back pressure 
is relatively high that the strength of the reflected shock is not ex- 
cessively strong (although it belongs to the strong shock solution) 
at- the triple poinr , and it soon degenerates .into a Mach wave as it 
extends from xhe triple point. It is then reasonable to assume that 
this variation may be expressed by (Fig. 33): 

= ' • ' “2 + ^"3 

. where cf^ is found from triple p'oint solution, 

is the Mach angle corresponding to , and 
Y and are the distances from the triple point normal 'to 

the flow as shorn in Fig. 34a and are nondimensionalized by 
using unit nozzle exit height Y^. 

Similarly, the- Mach shock profile may be approximaxed by 

where Y^ is the unknown Mach stem height and is -nondimensionalized 
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by using unit nozzle exit heighr, and 
0^^ is tile shock wave angle found from triple point solution. 
Since the airhient pressure is relatively high such that the 
central core f 3.0^7 will not be "choked" if it expands isentropicallj’" 
to the ambient pressure, it is thus possible to assume that both the 


upper and lower streams shall expand and tend tovrard the final asymp- 
totic rotational states corresponding to the uniform ambient pressure. 

A final balance of the flow momentum in the main flow direction associ- 
ated, with the control voltime as shown in Fig. 34a would allow a unique 
determination of the Mach stem height. 

From the principle of conservation of mass along stream t\ibe , one 
obtains . ' ” • . 


“2 “3a “ ''sa 


(54) 


and 


® “4a 


(56) 


where subscripts 3a and 4a refer to the asymptotic dovTnstream conditions 

at ambient pressure for regions and region 4, respectively. 

Using the maximum velocity V , local stagnation conditions, and- 

unit nozzle exit height Y as the nondimensionalizing parameters, one ob 

© , 

tains from Eqs. (54) and (55), 


dY 

•dY 




1/Y 
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u^ 
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'and 


dY 

dY 


4a 


(P.,/P .,)^^^ (uyv ) ^ol 

• 1 ol 1 m 


( 57 ) 
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where ^o4^^ol assmied shock profiles of 

Eqs. (52) and (53), and 



with similar expressions for u, /V , u^/V , and- u^/V as functions of 
^ 4a m’ r m’ 2 m 

pressure ratios . , - ■ 

The horizontal integral momentum balance for' the control volume 
shown in Fig. 33 gives 




'One now has three equations, Eqs. (56), (57), and (50), for the three 
mknovm cuantities Y , , and H, while Eqs. (52) and (53), in coh- 

junction with the shock relations', provide The resulxs of 
^o4^^ol corresponding stream tubes. 

Since and 6^ are given, with any initial estimate of one 
may find H from geometric relation 

2i' 

H_ = [1 - Y^ - tan 5_] cos 6. (60) 

2 o e 2 2 

where X is given in Eq. (50) (with Y =1 unity). Y and are 
e e , ■ o 2 



43 


subsequently divided into small divisions with each division represent- 
ing the height of a stream tube. Equations (56) and (57) can be inte- 
grated’ numerically using a trapezoidal rule to. find the unknowns 
and H and, simultaneously, the values of two unknoT-rri- integrals ap- 
pearing in momentum balance equation, Eq. (59). These values are then 

*■ t 

substituted into rhe momentum equation and the correct Mach stem height . 

would allow the integral momentum balance to be satisfied. This 
must be done in an iterative manner and. is equivalent to finding the so- 
lution to a nonlinear algebraic equation. The method of false position 
[29] works extremely well for the present problem. Upon representing 

the left-hand side of Eq, (59) as the residue function W(Y^), one starts 

1 - 2 

from two initial estimates of Y (identified as Y and Y ) such that 

o . o o 

1 V 2 

W( Y^) and U( Y^).have opposite signs. The next improved approximation 
will be - ‘ 


_ h W(^Y ) - h \-I(h ) 

3.^ _ o o ^ o_ 

W(^Y ) - W(^Y ) 
o o 


(61) 


3 1_ . 2 3 

If vK Y ) • vr( Y ) < 0 , one replaces Y by Y ; othen'fise one renlaces 
o o’. o o’ 

13 . ( 

Y^ by Y^. Using Eq. (61) for the next ’ approximation and repeating the 

calculation in this fashion until the criterion of accuracy is satisfied, 

that is; either! ^ I < e or j ^Y - ^Y I < e, where e is an arbitrarily 

small positive nturber, one then stops the iteration. 

A flov7 case of IL = 1.92, P ,/P = 3.67 (which is equivalent to 

1 ,ol a ■ 

^2 = 12.2071 degrees) is taken as an example. The two initial estimates 

1 2 . . “ 

are Y^ = 0.999 and Y^ = 0.001 vrhich represent two widely different Mach 

stem heights , and the final results are reached vjith three additional 
iterations. The actual result of all iterations based on the double 
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precision mode fr.om the I3H system 360/75 digital computer using the 
foregoing approximate analysis is presented' in Table 14 for various 
numbers of divisions N . vrith which H. and Y are divided. It agrees 
very \7ell with the picture given by A. Perri [43] from schlieren 
photograph technique- Another example of calculation- for = 1.92, 
3.54 (corresponding to - 12.9309 degrees) is presented in 
Table 15 which yields relatively lov/er Hach stem height in ccir.parison 
with A. Ferri's experimental data of Y^' = 0.4. The corresponding ve- 
locity profiles in the asymptotic dot-msti’eam flow are shown in Fig. 34 
for a variety of pressure ratios at = 1.92. 

Extensive calculations for two-dimensional overexpanded nozzle 
flows have been performed by this approximate analysis. The results 
of the height and the position of the Mach stem are shovm 5.n Figs. 25 


and 26. In general, rriis approximate method gives reasonabla results 


for flovf cases ’with subsonic state 3 at- the triple point. They also 


compare favorably V7-ith those from the derailed flo’.s' field analysis -irhen 
state 3 at the triple point is superson-ic. The cases corresponding to 


the tv?o extreme solutions of the triple point are also calculated and 


Identifie'd as before as 0„ = 90 degrees and a = 9G degress, it is in- 

’ O * V 

teresting to note that, at the extreme condition of a = 90 degrees, 

o 

Y /Y =1 and 'X /Y = 0 for high Hach number flow (M. > 2 . 5-) ; Y < Y 
oe oe- . l— os 

for Hach nunbers in the range of M, == i.5 - 2.5; and Y /Y - 0 for 

Which was predicted before in the triple point investigation. 
At the other extreme of = 80 degrees, exhibits a; maxi'rr.'um value 

in the neighborhood of M.j - 2,75, and bas a minimum value in the- 

vicinity of - 2.2. It should be noted that u'ithin the region .u'here 
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= 90 degrees, the flow is within the regular reflection flovr regime. 

Furthermore, because of the underestimation ' of the shock losses in this 

approximate analysis wifhin this flow regime, finite Mach stem heights 

are obtained from these calculations. Also shown in the same figures 

is the line of = 1.0 indicating that the stare at the triple point 

assumes the sonic state when it expands j.nto the ambient pressure. 

The analysis shovm in this section provides an approximate' method 

for estimating the Mach stem height in the region where M can be super- 
, ^ ' o 

sonic or subsonic at the triple point. Although the detailed flow field 
cannot be described, this method does give useful information in the 
area where the method of integral relations fails to yield any detailed 
information. The method is simple and- does not need extensive calcu- 
lations. Obviously, the extension of this approximate method to the 
'axisymmetric overexpanded nozzle flows is straightforward. 
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6. MACH DISC FROM lliBEREXPAi.'DSD AXISY^51-ETRIC NOZZLE FLOW 

6.1 BASIC DESCRIPTION OF THE FLOW FIELD 

For a uniform gas jet exhausting from an axisymmetric nozzle to 
a still medium at lower constant pressure, xhe flovr will first follow 
a tv70-diinensional Frandtl-Meyer expansion locally at the lip of xhe . 
nozzle. These expansion waves will eventually he reflected as com- 
pression waves from the constant pressure jet boundary.- When the 
ratio of exit to ambient pressure is high, the compression vmves of 
the same family vfill \mavoidably intersect each other .at a place where 
shock should originate to account for the coalescence of these vfavss. 

As the expansion and compression ^.-aves interact continuously downstream, 
the shock becomes stronger because of the accumulating effect of the 
coalescence of these reflected waves. This iinbedced shock wave could 
be reflected from the centerline of symmetry as shoi-m by Pack Cl2] for 
a two-dimensional gas jet. However, if the surrounding pressure is • 
low enough, the snowballing compressive effect will result in an im- 
bedded shock wave at downrstream field so strong and curved that regu- 
lar reflection from the centerline of symmetry is impossible . The con- 
figuration associaxed with Mach disc xhen comes into play. 

Figure 3 shows a typical Mach configuration from axisymmexric under- 
expanded nozzle flow. The previous discussions for the t-.-ro-dimensicnai 
overexpanded nozzle flow are applicable xo this axisymmetric under- 
expanded problem bux with several basic differences. First of ail, the 
incident shock is no longer straight, and the flow upstream of the re- 
flected shock (region 2), is highly nonuniform.. For simplicixy, the 
central core flow' behind Mach disc can be approximated by a one-dimens icr.al 
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analysis which has been used successfully in [15] for overe:<panded 
cases. Since the incident flow in front of the Mach disc is also non- 
uniform and not horizontal locally at triple point, a rotation of co- 
ordinate as shorn in Fig. 4 for triple-point solution is necessary. 
Furthermore , the rotationality is no longer negligible because of the 
■ rapidly changing shock strength behind the incident and the reflected 
shock waves. The method of characteristics will be used exclusively 
for the calculation of flow structure of the supersonic axisymmetric 
gas jet, while, one-dimensional isentropic flow is applied to the 
central core flov? immediately after the triple-point shock system is 
inserted into the flow field. 

7^6.2 THE METHOD OF CHARACTERISTICS 

. In axis\^mmetric, supersonic flow field the fundamental equations, 
'-Eqs. (14), (15), and (16), in cylindrical coordinates can be solved by 
“-following characteristic curves in physical and hodographical planes 
s-i-multaneously [31,32,393. The characteristic equations are 


dr 

dx 

= tan 

(0 + a) 


■ (62) 
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The added last term expresses the difference betv/een axisymmetric and 

two-dimensional problems. Application of these equations to the field, 
free jet, and axis point procedures is well known and is given in AP- 
PENDIX C. 

During the process of computation , the characteristic curves of the 
same family in physical plane are foiuid to intersect each other. It 
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Indicates the failure of the continuous, single-valued solution, and a 
shock discontinuity should accordingly be introduced. into the f low • 
field'. This is where the imbedded shock v?ave originates . The reflected 
wave resulting from the intersection of the compression wave and shock 
wave should not be neglected since it will reflect from the free jet 
boundary and has the augmenting effect on the dovmstream shock strength. 
In this Mach disc problem, the accurate calculation of incident shock 
strength is^ fo’und to be critical since all the downstream flow field, 
including the triple point solution, depends on it. Therefore, the re- 
gion where compression wave and shock wave of the same family intersect 
to produce a stronger shock wave should be examined carefully and is of 
paramount importance. Moe and Troesch [41] are correct in stating that 
^^^the calculation of the imbedded shock wave is the most complicated part 
•of the numerical method. The simple fold-back method [16] or other 
simplified methods, though determining the shock position properly, 
will not give the appropriate shock wave strength . 

As shown in Fig. 35a, I is the characteristic of family 1 and II 

is the characteristic of family 2. The properties at E, A, and C on 

characteristic curve II , and at F on II are knovm from previous cal- 

culations. The flow variables at B are found from known properties at 

A and r by following characteristic curves AB of I and FB of II . 

m n 

Similar calculations can be formed to. find point D with DC of and 

BD of II , However, if the x-coordinate of D is detected to be at the 
upstream part of 3 on the characteristic it indicates waves AB and 

CD have intersected each other somewhere at G as shovm in Fig. 35a. 

Whenever waves of the same, family intersect each other , extreme 
care should be taken to determine the flow properties at the intersection 
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point. As shovm in Fig. 35b, thei’e appears a slipline representing the 
difference in entropy beti^een and and. a reflected rave RG ac- 
counting for the pressure difference at G„ and’ G, , . The coordinates at 

3 4 ‘ 

G are deternined froti geoTr-etric relations of the intersection of two 
straight lines, AB and CD. The flow properties (H and G) at are 
established from a characteristic calculation v;ith point A and a point 
between E and F. The flow Each nuinbers at between vraves AG and CG 
and at G^ behind wave CD are obtained from linear interpolation of the 
flow variables at A and B and at C and D, respectively’'. The oblique 
shock \;ave relation is used for both compression wave and shock v?ave. 
Since compressive vjaves converge i one has 



here, 6 is the deflection angle found from and across the shock. 
Since reflected wave RG is usually very weak, one may treat it as 
an isentropic wave. The Prandtl-Heyer function is then employed to 
find the flow angle at G , from nroperries at G_. The flow variables 
at Gj^ behind the merged wave GP is determined from oblique shock re- 
lation . An iteration’ procedure is then necessary to determine the Mach 
number and flow angle at G, and G,, , from the condition of equalitv of 
flow angle and pressure across the slipline GS, which is in a very simi- 
lar situation as- that encountered earlier in triple point solution. 

The shock is assumed locally straight Vfith constant strength extending 
from point G up to point P where it intersects with characteristic 

curve II and replaces both points D and B in the flow field. It is to 
n ■ ■ 

be noted thax the variation of this shock strength due to wave coalescenc 
is accounted for in this manner and the existence of the interface 
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.63 is completely ignored afterward. 

This is a systematic manner of treating the wave intersection 
problem since the resulting shock position and strength must be 
registered and recorded at all times. The flow angles at and 
cannot be determined from linear interpolation as was done for their 
Mach numbers '.since the local shock relation might be violated al- 
though the difference between these calculations is very minute in 
all of the flow cases -with small characteristic grid. Furthermore, 
the reflected wave RG inserted into the flow field has the effect of 
reducing the grid size and increasing the accuracy of computation. 

6.3 THE CALCULATION OF INTERACTING FLOW FIELD 

Based on an idea similar to that in the ti-fo- dimensional over- 
expanded nozzle flow, the- position of Mach disc is unknovm and must 
be solved by an iterative fashion . The correct location and the radius 
'of the Mach disc are determined from the considerations that as a re- 
sult of the interaction be-tween the waves and the_ flow fields, the cor- 
responding uniform one-dimensional analysis will give a sonic condition 
a throat dovmstream of the triple point. 

The imbedded shock strength from the method described above is 
checked at every step of calculation toward downstream to see whether 
or not it satisfies the condition of triple point solution. When this 
condition is met, the central core flov/ is then approximated by a uni- 
form one-dimensional trea'tment and the flcr.r field Including the reflected 
shock v?ave s'trength in the super>sonic upper part of the slipline is 
solved by the method of ‘characteristics . Again,' the Mach disc stand-off 
distance is assumed to be zero as has been done in the •two-dim.ensional 


case. 
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• The very first step along the slipline is taken' to be very small 
and to be a straight line having the slope from the triple point so- 
lution* As shown in Fig. 36, the first segment of the slipline TS has 
the slope tan 6 where 6 is found from rotation of coordinates 
(0=6^- 0^ given in Fig.' 4). From the condition of equal pressure along 
the slipline, the flow Mach numbers and %1 at the end of the first 
step can be solved from the following isentropic and one-dimensional 
equations : j . . 
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and 

^ = r + (x, - X ) tan 0 
1 o 1 o 

. The characteristic net layout for the subsequent step-by-step cal- 
culation is shovm in Fig. 37. The flow properties are knovm at T and 
from previous results. The flow variables at where characteristics 
of family II and the reflected shock intersect are found from solving the 
characteristic e'quations for family II and the oblique shock ^rave 
equations simultaneously as presented in APPENDIX C. This must be done 
in an iterative manner. • One added unknovm during the iteration procedure 
is the flov; condition in front of the reflected shock wave whi ch is 
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obtained through an interpolation procedure. -As shovm in Fig. 38, if 
lies inside the characteristic grid ABCD, then the Mach nuniber and 
flow angle in front of the reflected shock can- be approximated from 
linear interpolation ; 


M- 
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(67) 


Near the free jet boiondary the grid can be triangular Instead of 
quadrangular, and a corresponding modified form consisting of only- 
three characteristic points should be used. 

The need for the properties in determining the flov? variables up- 
stream of the reflected shock wave' indicates one of the reasons that 
the characteristic nen must be stored and registered. The method 
shown above for non-uniform approaching flow is the extension of the 
well kno™ technique for determining the shock wave location shov.ti in ■ 
Fig. 17.21 [31] for uniform free stream conditions. 

■ Figure 39 illustrates how the location where reflected shock inter- 
acts with the free jet boundary is determined. VJhen the coordinates of 
the point on the reflected shock are found to lie outside of the free 
jet boundary, the point of inuersection between the shock and the free 
jet is found from the intersection of two straight lines, AB and CD, as 
shown in Fig. 39. The flow angle in front of s'nock \rave is found from 
linear interpolation between points A end B , and the Mach number is that 
of the free jet boundary. The oblique shock wave relations are then 
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employed to obtain the flovf properties behind the reflected shock wave 
subsequently followed by a P-M expansion. 

Similar to that of a two-dimensional problem, the slip line point 
procedures are governed by the following equations : 
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where Q(M) is given in Eq. (45). The simultaneous solution or these 
equations will give the flow variables across the slipline. 


6.4 RESULTS AND DISCUSSIONS 


The analysis presented above for the underexpanded axisymmetric 
nozzle flow is programmed successfully on an IBM system 360/75 digital 
computer using double precision mode, A single successful run of The 
Fortran IV program takes less than four minutes compuTation time. For 
the convenience of reference, Fig. 40 illustra’ces the free jet kach 



numbers and press\ire ratios . Also shovra in the same figure is the cor- 
responding ratio of ambient pressure to nozzle stagnation pressure at 
different nozzle exit to ambient pressure ratios. 

Figure 41a illustrates the results of iteration for determining 

the Mach disc location for the case of M = 2.0, ? /? = 6.39 (? /P = 

e e a a oe 

0.02) disregarding the effect of vorticity for = 20 where N denotes 
the number of uniformly spaced vraves initiaxed from the initial Mach lin 
Tne early insertion of the Mach disc, for example at x = 11.86, will 
cause the slipline bending continuously upward yielding only subsonic 
states in the downstream flow field. On the other hand, the insertion 
of the Mach disc at x = 12.025 will cause the central core flow reach- 
ing a, sonic condition where the area is not a minimum. The correct lo- 
cation of the’ Mach disc for this flow situation stands at x =; 12.00. 

-It is obvious that the flow field is very sensitive to the location of 
the Mach disc; a slight change in the location of the Mach disc will 
result in a tremendously differenx flow pattern downstream. 

It has been mentioned previously .that because of the enormous 
difference in stagnation pressure behind the curved imbedded incidenx 
and the reflected shock waves , the roxationality cannot be neglecxed 
in this study. Indeed, inclusion of the vorticity into the consider- 
ation results in a sizable shift in the location of the Mach .disc. 
Figure 41b gives a comparison of the numerical results of the flow's 

with and w'ithout vorticity for the case of 1-5 = 2.0, P /P =6.39 

e e a 

(P /P =0,02) for N = 20. These results clearly demonstrate the in- 
a oe 

adequacjr of the treatment ignoring the vorticity v;hich gives a much 
smaller Mach disc standing farther av;ay from nozzle exit. In all of 
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the follo^fing results, rotationality has been included in the consider- 
ations . 

These improved calculations with vorticity, however, exhibit a. 
somewhat unexpected phenomenon. In the determination of the -sonic 
state in central core, xhe flow angle does not vanish and thus the 
area is not -a minimum. Moving the Mach disc toward downstream will 
resulr in a steeper flow angle at sonic condition. On the other hand, 
no triple point solution is possible for early occurrence of the Mach 
disc. Presumably, this is because of small inaccuracies in the up- 
stream numerical calculation. The reason is tvrofold: First, the 

linear interpolation used here in the characteristic grid is not accurate 
enough to handle the rapid change in shock strength near the triple 
point where the flov? variables depend sensitively on the incident shock 
strength. Secondly, .the change in stagnation properties through the 
variation of the entropy gradient term is based on the - interpolating 
.properties from known states, and is, therefore, dependent on the lo- 
cation of streamline passing through the unknovm state (see APPENDIX C). 
IJhen the stream.line location and direction at the unknovm state have 
different values from those found from linear interpolation, the entropy 
gradient term x^rill not be calculated correctly. This influence is more 
pronounced at or near the occurrence of the triple point v/here the de- 
flection angle of the shock changes rapidly. Taking these circumstances 
into consideration, the results presented for flow with vorticity will 
be considered as the upper limit of the location of the Mach disc. 

A typical characteristic vrave pattern in physical scale from the 
actual computation of = 1.5, f 5.0, N = 30 is given in Fig. 42, 
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The location, and strength of the imbedded incident shock resulting from 
the coalescence of the compression waves are computed automatically by 
digital computer from the method described before whenever waves inter- 
section occurs , 

In all the calculations of this study for M =1.5, 2.0, 2.5, and 

e' 

3.0, the sonic state of the central core flct-7 is reached before the ap- 
pearance of the P-M fan which results from the intersection of the re- 
flected shock with the fr-ee jet boundary. The computation is stopped 
at the place where the sonic state of the lower side of slipline is de- 
termined. This explains the unfinished calculation of the reflected 
shock and its dcnmstream flow field which should be in a supersonic 
state and does not have any influence on the already established Mach 
disc and -upstream flow field -although the calculation .of v^hole down-- 
stream field including P-M fan is well programmed and is incorporated 
in the same Fortran code . 

Some of the free jet boundaries and shock wave patterns for M^ = 

1.5 and 2.5 are shown in Fig. 43 and Fig. 44. It is evident that higher 
pressure ratios (P /? ) for a given nozzle Mach number would I’esult in 
larger Mach discs at farther downstream locations. At the same pressure 
ratio, however, a reduced nozzle Mach number brings Mach disc closer to 
the nozzle exit plane. 

Figures 45 and 46 show the distributions of various flow .properties 

for N = 30 at M = 2.0, P /P = 3,195 (P /P = 0.04), and at M = 3.0, 
e ’ e a a oe ’ e , 

P /P = 5.0 (P /P = 0.005445), respectively. All the sxagnatior. nres- 
6 d 3 O0 

sures are normalized by nozzle stagnation pressure. The stagnatic- 
pressure behind the incident shock drops rapidly near the triple p 
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which is a consequence of the accumulated compressive effect on the gas 
jet. Immediately behind the reflected shock, however, the stagnation 
pressure increases and deviates away from the constant stagnation pres- 
sure on the upper side of slipline. All of these signify that the 
rotationality certainly needs to be included in the considerations. 

The differences in Mach numbers and in flow angles betT'fsen the 
states in front of the incident shock and the centerline of symmetry 
decrease along the flovf direction and are very close to zero in front 
of the Mach disc. This demonstrates that the flov; in the very front of 
Mach disc is rather uniform for the cases shown in Figs. 45 and 46. On 
the contrary, the differences in Mach numbers and in flow angles in. 
front of and behind the incident shock are obviously increasing toward 
downstream illustrating the growth of shock strength -and, eventually, 
reachxng the state where the trxple poxnt solution is able to be satis- 
fied. . , 

The -trace of the imbedded shock appears when the reflected compres- 
sive wave from the free jet boundary first intersects vfith the last P-M 
characteristic wave from the lip of the nozzle. In this numerical study, 
the position VJhere discernible shock wave originates is given in Fig. 47 
when P-M fan at the lip of the nozzle has been divided into approxim.ately 
one-degree increments. The location where the imbedded shock wave has 
the maximum radius is shown in Fig. 48. Figure 49 gives the physical 
coordinates where the free jet boundary has the maximum height. Shown 

in the same figure are some calculated data at P /P = 10 from £16] us- 

6 a. 

ing the fold-back me-thod which agrees perfec-nly with the results of the 
present investigation. At the same pressure ratio higher noz- 

zle Mach numbers alvrays produce wider gas jet profiles with maximum 
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width occurring farther downstream. 

The Mach disc location and radius are presented in Fig. 50. Some 
experimental data from Love and Grigsby [16] are also shown in the same 
figure for comparison. Since the characteristic grid size is finite in 
the numerical calculations , it is , therefore , apparent that the actual 
‘ Mach disc will occur earlier because of the faster growing shock 
strength associated with the intersectipn of infinite number of waves. 

For the same reason, the Mach disc of the numerical calculation appears 
_ shorter, than the actual one does. For higher nozzle Mach number and ■ 
^■-far away from the initial characteristic line, the characteristic grid 
..size, becomes quite large which further explains the overestimation in 
„.the' x-coordinate of the Mach disc location. On the other hand, at 
. lower pressure ratios, the Mach disc is very small and the characteristic 
net near centerline of symmetry is- usually not 'fine enough to produce 
--good results near the triple point and tends to underestimate the radius 
. of the Mach disc. Both effects ‘are more pronounced at 'higher nozzle 
-r-Mach numbers 

Table 16 illustrates the effect of changing "the number of initial 
..grid divisions' (N) on the calculated results. The larger value of N cor- 
■- responds to finer initial grid size. As indicated in the table, the use 

- of larger N does not necessarily give better results compared with experi 
. mental data. The -finer the initial grid, the more- iterations are needed, 

- and the accumulated numerical errors might aggravate the final ansv-er. 

All of the calculated results show that N = 25 or 30 is a fairly good 
choice for the nuni)er of .divisions on the -initial characteristic line. 

Another manner of getting more accurate results is the insertion of 
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extra field points into rhe flow field whenever- the grid size is de- 
tected to be larger than a predetermined value . This will create 
smaller characteristic grids in dovmstream field without imposing a 
large value of N on the initial characteristic ].ine and will probably 
eliminate unnecessary growth of numerical errors. However, this will 
somewhat complicate the programming work especially irhen all the wave 
coalescences are to be considered. 

One of the salient features in most of the calculations is that 
the slipline tends to go upward first in the vicinity of triple point 
and then bends dovmward decreasing gradually to reach a. minimum height . 
at the throat (for example, see Fig. 42). Figure 51 shows the approxi- 
mate coordinates of the throat as found from numerical results -when the 
center cere flow behind the Mach disc is close to sonic state. 

Finally, Fig. 52 reveals the information necessary for finding the 
properties at the triple point- The resulxs are so given that the triple 
point solutions of Figs. 5 through 9 -(or Tables 1 through 13) can be 
found easily from the corresponding and 6^ given in Fig. 52. From 
the computed data, it is noted that the Mach disc occurs in the very 
neighborhood of the limiting case of curve 6^^ in Fig. 10a. Unlike the 
two-dimensional overexpanded flovf, the incident flow angle 6^ (see 
Fig. 3) is not zero at. the Triple point in these undefexpanded gas jets, 
and the interpretation of the angular relationship associated with the 
triple point solution must take this inclination, 0., into proper account. 
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■ 7. CONCLUDING REMARKS 

Based on -the numerical results obtained through extensive calcu- 
lations on problems with Mach reflection or Mach disc; it may be con- 
cluded that: 

1. The one-strip method of integral relations and the method of 
characteristics , together with the one-dimensional isentropic 
nozzle analysis for the downstream part of the central core, 

, are usefiil to give a more detailed description of the Mach re- 
flection pattern from the two-dimensional overexpanded nozzle 

J ’ 

flows with weak reflected shock at the triple, point . Their 
application to strong reflected shock is hot successful. 

2. The approximate method to estimate the Mach stem height for 
•the tiTO- dimensional overexpanded flow problems , which was 
originally developed for flow cases vrlth strong reflected 
shocks, would give meaningful results even if the reflected 
shock at the triple point were weak. 

3. For axially symmetric gas flow -with Mach disc, the vorticity 
must be taken into account . The wave interactions , including 
their coalescence into shock waves v^ithin the jet flow field, 
are extremely complicated; accurate calculations seem to be 
difficult, if not impossible, to achieve. 

4. Although the viscous effects tend to modify and influence the 
ultimate flow patterns, it is believed that the inviscid in- 

• ■ vestigation of these problems is already yielding fruitful re- 

sults . 



FIGURES AND TABLES 




Figiire 2a Regular Reflection produced from Two- 
Dimen'sional Overexpanded Nozzle Flow 



Figure 2i> Mach Reflection produced from Two- 

Dimensional Overexpanded Nozzle Flow 







Figure 3 Mach Disc produced fr'om Axi symmetric Underexpanded Nozzle Flow 
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Figure 4a Limiting Case of Triple Shocks 
Intersection ■ = 90°) 



Figure 4b Limiting Case of Triple Shocks 
Intersection = 90°) 















REGULAR REFLECTION POSSIBLE 



Figure 10a, Regions of Regular and Mach Reflection, vs. 
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Figure 10b Jet Deflection Angle S vs. Pressure Ratio P„/P 
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Figure 11 Oblique Reflection of Shock Wave from the Central 
Plug at a Nozzle Exit 


Y Y 



Figure 12 Detailed Figure for Mach Reflection from Two-Dimensional Overexpanded Nozzle Flow 



















Figure 19a Illustration of Applicability of the One 
Strip Method of Integral Relations for 
Weak Attached Shock Wave 
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Figure 19b Illustration of Inapplicability of the 
One-Strip Mexhod of Integral Relations 
for Strong Attached Shock Vfave 












Figure 23 In-Scalc Plot of Two-Dimensional Mach Reflection Pattern (M, = 3.0, 6 
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Figure 25 Mach Stem Height Y /Y vs. 6„ for Two-Dimensional 

o e 2 

Overexpanded Nozzle Flow 
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Figure 33 Control Volume and Coordinates used in the Approximate 
Method for Estimating Mach Stem Height and Location 
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Figure 34 Velocity Profiles in the Asymptotic Down- 
stream Flow from the Approximate Method for 
Different Pressure Ratios at M, = 1,92 






FigTore 36 First Step along the Slipline for 
Axisymmetric Under expanded Flow 
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Figure 37 Characteristic Net Layout behind the 
Reflected Shock Wave 
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Figure 38 Determination of the Flow Properties in 
Front of the Reflected Shock 



Figure 39 Determination of Flow Properties and 
Coordinates at the Intersection Point 
of Reflected Shock and Free Jet Boundary 
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Figure 40 Ambient-to-Nozzle Stagnation Pressxire 

Ratio, P /P , Free Jet Mach I'lumber, and 
’ a oe ’ r’ 

Initial Free Jet Inclination Angle at Nozzle 
Lip, 0^, vs. Nozzle Exit-to-Ambient Pressure 
Ratio, P /P_ 







Ila DcpcndeiiC'j of Hicj IJonic Gtato on die Iioccition of Mach Disc v/ithout ];Til..'>ony Gxv.dicnt 
at H = 2.0, P /r “ 0,30, K = 20 
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ROTATIONAL ANALYSIS 

ANALYSIS WITHOUT 

VORTICITY 


2 3 4 5 6 7 8 9 10 i I 12 13 14 15 


Figure 41b Comparison of Analysis with and without . 

Vorticity atM =2.0,P/P =6.39, N= 20 
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Figure 42 


In- Scale Charahteristic :Wava Pattern in 
Physical. Scale at = 1.5, ^ 5.0, N = 30 
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Figure 43 Free Jet Boundaries and Shock Wave Patterns 

at Different Pressure Ratios for M =1.5 
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Figure 48 
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Coordinates of the Maximum Radius of the Imbedded 
Shock Wave, X, (at r ) and r„,, , vs. P /P 
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1 . 7 08 

4.043 

2.802 


17.000 

0.34 7 

1.2 79 

0.836 

0.641 

‘r . '1 3 9 

12. ‘‘cl 

'^0,5:. 9 

79.590 

81.103 

2.398 

1.681 

,031 

2.3-07 


1 / .250 

0.352 

1.297 

0.83 7 

0 . t " 7 

V . 0 1 0 

13.210 

50.990 

80.21-'. 

00.536 

2.-3 0 

1 . 6 3 

<•-.018 

2.812 


17.500 

0.356 

1.255 

0.839 

0 . o 5 - 

3.63 0 

13.870 

51. 431 

80.895 

79.937 

2. 44 2 

1.6 26 

- . 003 

2.813 


1 /, 750 

0.361 

i. # 2 'V 3 

0 . 34 1 

0 . 77 1 

3.208 

I'f. ;■ 4 2 

51.8 8<' 

81.642 

79.300 

?.6 9r,. 

1 .5 97 

3.98c 

Z-SZ-”' 


1 B . 0 1 0 

0 e 3 O 7' 

1.230 

0 ♦ 

0. c 70 

2,7 75 

15.225 

52.351 

82.46 7 

7 0 . t 1 8 

2. *=29 

1.56 8 

3.96 7 

2.831 


10.250 

0“. 3 71 

1.21 7 

0 ,8'^. 7 

0.6 80 

2.33 0 

15.920 

=2.831 

83.380 

/ 7 • b O 

2 . f'64 

1 .6 38 

2.945 

2.840 


18.500 

0. j76 

l,2 0-> 

0.85 1 

0 .691 

1.8 7‘ 

16.626 

5 3. 32 8 

S-. .400 

77.085 

2.6 00 

1,607 

3 . -U 9 

2 . 8 •• 0 


10. 750 

0. 382 

1.190 

0. 65 5- 

0. 7 0'- 

1.4 0-; 

i 7 . 3 4^ 

5 3 . 8'~2 

8: .8-i^9 

76.208 

2. 637 

1.47 5 

3. -19 0 

2 . 8 6 1 


19.482 

0.3 98 

1 . 1 V7 

0.877 

0.75 8 

0.000 

19.482 

55.479 

90.000 

72 .94 9 

2.753 

1.367 

3 . 75 4 

2.'509 
















H 
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Trip le-Point Solutions , M, = 2i0 


^n = 2»oon 




OlILTA/V 

'o'. 0 


l.YS'O 
2,268 
2r7ii ■ 
3.310 


‘i2.3r2 
J, 3. 01)3 
43'r3'‘-;=i' 
43, 710 
T4 . 021 
4-'-.33i 


378 

ooy 


ill 

. 3G3 


S1G‘'!A3 

70, V 73 

71. 'L97 


71 . 

71.j<-2 7 

TF. Fm" 

71.7 0? 


7 1 , 3 " 8 
72.0^2 


72,232 
72.‘,' 3 7 


310!^ A 4 
90.000 
89.232 


UH ,9‘Mi 
3 8,63 8 
88.322 
86,004 


3 7.88 
87,33^ 


b7.02i 
86 « b 1 


P2/P1 

1.9‘'9 

2,011 


2.082 
2 . 08 8 
2.1 L-' 


2, UO 
2.167 


2- 196 
2 , 222 


P3/P2 
2 . 3 yf 
2.2 38 


2 

2.181 

2.15^. 

2.126 


2.099 

2.072 


2.0-6 
2 . 0 1 0 


P- /P 1 


4.6 00 
./-99 


(P4/PI ) 
2.980 
2.9 80 



2 . 986 
2.986 


6. 006 

"46 . 02 7 

72. 360 

iTo ,335 

2.250 

1 . J 9 1 

'+-■'-81 

2.90 8 

6 . 5 c 3 

•O': .376 

72.900 

8 ’" .982 

2.2 79 

1 .o4c: 

4.777 

2,9 89 

7, 126, 

46-731 

73.159 

>J-.62i 

2 .'3 08 

1-938 

7f . 7 3 

2.991 

7- 695 

7- 7,. 090 

73-639 

05. 2 "^2 

2.357 

1.912 

7- . 4 8 

2 . 993 

8 .'2 72" 

"4 7 .'Vi 5"" 

73, 7 ^0“ 

“84.872 

2.36 6 

1.086 

4,4.6 3 

2.995 

8.855 

4 7. 82 7 

76, 065 

84.483 

2.397 

1.8 50 

4.457 

2.993 



13, 173 
U. 830 

'i4.r97 

13. 1J7 
"ll". 8 70 

1 . 6 7 r- 

“l 7.29 7 
18.0 32, 

ioFfoi 

19. 643 
21,329 


502626 77. 
5l.0.t0 77. 
"5 1 . Fot" ' 73 ; 
51.9^4_79, 
'5 2.4’F5 79 ■, 
5 2 . ^2 1 30, 
■s 7F2 3 (il, 
5 3.90 32. 

84. 

Oi.Of-O B5. 
' ['ll. 0 79 90. 


.35 3 
',325 

r.2 6 7" 
7o-vr 

i . 3 6 5 


'.6 30 
.825 
.~9l-7 
.9 VI 
,93’J 


2.622 

2 . 6 92 
2, 723_ 

2 . 804 


2 . o => 3 
2 • 8G^ 
~2."92 5 
2 ,9-'0 

■ 3 . 07 / 


1.976 
1.60 
I ,"/22‘ 

J_._0 9^ 
1 , 3 7. 6 
1.337 


l.-JO.s 

L*F IF 

i . i 4 ic. 

\ .•' ■10 

1.317 


4 . 3 O'" 
4,331 
4 73 6 7 
4.3 6 0 
^.331 
A .310 


4.2 86 
■4. 250 
4 .'2 2 > 
•..185 
- . 0 '■ 1 


3.019 


3.02''- 
3. 02 9 

3.03r 

_3. 0<''2 
3 V 0 PO 

3. 0=^9 


3.069 
3 . 0 0 l_ 
~3 7o9^ 
„3 , 1 1 2_ 
3.169 
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Table 5 Triple-Poi nt S olutions > = 

. I .. 1 I ■ ( I I ' . ! ■ 1 ' ■ > H ' I d’ . . ■ , 

Ml = 2.200 


TFIPLC POINT SOLUTION 


GAMMA = 


DELTA 2 P2 /P 01 
~TS,r2u"TT7n6 L 
15.300 0.217 1 
~l3T6'0-0“)".22rT 
_15.«->00 0.22^ I 

~rcT2 OTTTZZiTT 

16.500 0.231 1 

■"■r6';aoo"o,235'”i 

_J.7.ino 0.23» 1 
i"7l‘^,-0?r 0.2^2“ i 
17.730 0.246 1 
1 ir.TiT) (TTj . 2 Z y r 
18.333 0.253 1 


M2 M3 
.fc'20 oToiff* 
.612 0.918 
ToO'O TT.'oVJ' 
.387 0.918 
372“3“."5'nr 
.5*32 0.91/ 
.5A:9”'0.'917' 
.536 0.916 
■.52 3'"0.9U'' 

.510 0.915 
.A yo'‘07‘91-r 


M4 

T3~."5T/' 
0.517 
ir. 547 

0.548 

0 . 5 ^‘9 
0 : 549" 
0.5 50 
0'.5 5'2" 
0.553 
(TT 5 9 
0.5 56 


DELTA 3 
"i 5'Vi2 8 
14.935 
"iV.'orro 
1‘i-. 339 
'i“cir4 
13.703 
13.3/6 
13.044 
'i2'. IQl 
12.362 
"1 Z'.O'l 1 
1 1.654 


_Dt;LTA4 

o;o" 
0.345 
~ 0 ; 9 5 if 
1.561 

2 . i Vt" 
__2_. / 97_ 

3 ". 4 24"’ 

, 4. 0 58_ 

—■4 . 94" 
_5 ._3 3 ^ 

5"1 9f39 
6 . 646 


5 I_OMA 2 
4 1 . 4'22 
41._p30 
"41'. 99‘5’ 

'n2". ?Ty 
*3.113 
■4 3 , ?/02" 
4,3.89_1_ 

44,236 
■ 44.686 
■"( 5 :'oT;2“ 
4 5.5 05 


S 1 r.MA3 

S IC.MA4 

P2/P1 

P3/D2 

P4/P1 

1 P4/P1) 

t;5. 75* 

90.000 

2.306 

2.37 7 

‘“.480 

3 . ■" * 8 

65.7 94 

89.823 

2. 3 2.5 

2, 38 7 


3.5 48 

66.875 

89.613 

2 . 3 6 1 

2.321 

6, ‘■•00 

■KWBHMI 

u5 .970. 

89.199 

2.397 

2.286 

"■.4 79 

3 9 

6 6 . 0 79 

88.383 

2.* 3* 

2, 2-^0 

5 , T 78 

3 , M 4 y 

66 .204 

88-562 

2 .471 

2.216 

5.4 76 

3.650 


v/r;34i 

66.502 
■6t'.6"73' 
56.871 
i 7.0 '8 5 
67.319 


¥b:2 3 8’ 
87i908 
8 7. 57 ‘4' 
jl/_^2 3 5 
8L.aa'9" 
86.536 


2.5/./ 
'2. 5 or 
2.626 
"2 . .6 6 6 
2. 706 


‘18'2 r,'5 5’l 


2. 1-.-8 
'2 .1 1 5" 
2. 082 
2 .'049 
2,017 


^.i_72_ 
7 0 
5 . 4 5 7 
5 - ' * 3 
5.45 9 


3 . 5 " 2_ 
3 .'5 53 

3.555 

3.^57 

3.550 


3 . ' 6 1 
3.‘=63 



2 I . 6 00 
21.900 

" 2 ' 2 '."200 ■ 


2 3. 700 


0.299 1.3 23 0.895 0.-.96 /.131 i*.-i'69 

0.304 _1 ,3 07 J').893J3.6 )2_ 6.o4_7 15_.253 
Tf. 30 9 1. 2 9'0” o'; 89 1 n ; 60'9 “ o" l7, ,f “i> .05 " 

0.^1-'- 1 .273 n. 889_0. = 17 5.(:27 16.37^ 

o"'.3i9”i .‘zsirnTuda alelt T2o'd7'~u 

0.325 1.2 i>3 0,887 0.d37 4.523 18.5/7 


0.350 1.220 0.886 0 
0.336 1.201 O. 3 O 0 O..t.o3 3.318 2'> 


'■0.5/2 
. 0'-;0 
f? 1 , r> '■ 0 

52.200 


52. 779 
':3. 3 80 


30 7 
r>6 3 


71 .838 
72^5 IJ^ 
■7 3.25‘. 
74.0 7.3 


74 . 9 9^ 
76.026 


77.19'- 

78.535 


8 1 .904 
81.368 
" 80 '.”8 no" 
80.193 


79-539 
7-8.82 9 


78.048 

77.177 


3.2 0? 
_3.2:>3 
3.‘'30'5" 
3.359 


3.^ 

3 . 7 1 


4/6 

5 . 3 c 8 

3 .6 05 

^45 

6.353 

3 .-'-13 

>"l'4 

5.334 

3762 r 

^03 

5 . 3 U-. 

3.^31 . 


5.294 

3.64?. 

: 18 

S' . 2 fc y 

. 6 " 


3.690 


24 . 000 

0. 3 '-2 

1.181 

0.88 7 

0.77 80 

2.669 

2 ! . 331 

5*^ ,36 6 

80.095 

75.188 

3.656 

1.411 

5.158 

3.713 


24.300 

0.3‘»8 

1.159 

0.890 

0.701 

1.98 7 

2 2 . 3 1 3 

5C. 092 

81.9 

76.036 

3. 723 

1.371 

: . in.’. 

3. 74 ;’ 


2-» .C'OO 

0.35 5 

1.137 

0 , 893 

0 . 7 2 7 

1.271 

2 3. 32 9 

56.883 

0^' .2 38 

73.447^ 

3.794 

1,326 

6,032 

'3 , 7 8 0 


2 5.119 

0 . 3 7> 0 

1 .09 J 

O.vl 7 

0.8 00 

0.000 

2 '7 1 1 9 

I- r, . - 1 

90.000 

70.174 

3*. 9?:- 

1 . 228 

4 . 8 3 '1 

3 . 8 8 8, 
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' Table 7 Triple-Point Solutions, M^j=_ 2.6 


TRIPLE POINT SOLUTION 

DELTA2 P2/P01 M2 MJ DELTAS DELTAA 

L8.i41 0.151 1-bO'J i-OyO 0.50A 10.141 0.0 

18.400 0.153 i,79o 1.079 0,504 17.890 0.510 

- 10.800 0.156 1.778 1.077 0.S04 17.499' I-'JOr 

19.200 0.159 1.758 1.075 0.505 17.103 2.097 

I9.60a-01 153'"l-739-l. 072-0, 505" 16;‘703 2. 897“ 

20.000 O.lo6 1.720 1.069 0.506 16.297 3.703 

20. ‘'-00 0;i69”i.700 l-066''0.507 1 5. 8 Bo'"' 6, ?1 4 

20.800 0. 173 1.681 l.Oc-2 0.508 15.469 5.331 

21,200 0.177 1.661 1.058 0.5 10 15 .0-<6 -‘o- 154 

21.600 0,180 1.641 1.054 0.512 14.616 6.984 

22-0a0''0;i34~i,621' "1.050-0.5 14 '14:178 7,622' 

22-400 0.188 1.600 1.045 0.516 13-733 8.667 

— 22-800 0.192 1.579 1.040 '0.519 13.278 - 9,522 

23.200 0.196 i-558 l.0'34 0.522 12.814 10.386 
23,600 0.200 1.537 1.029 0-525 12. 340 "11, 260 

24.000 0.204 1.516 1.023 0.529 11.854 12,146 

24;400"‘0;208“-i.4 94' 1. 0i6“'0;'533"i l':354""13. 04'6' 

=,.800 0.212 1.472 1.010 0.538 10.041 13.959 

25.200 0.217 1.449'1.003 0.543 1 0 - 3 1 2 ' 14 . 8 88 

25.600 0.221 1.426 0.996 0.549 9.765.15.835 

-”■"26,000 0.22o 1.402 0.989 0.555 ' 9.198 16. 802' 

26.400 0.230 1,378 0,981 0-562 __8.600 17.792 

2o. 800"'0.235“1.3'54 ■0.y74"‘0.5 7i ■“ 7 ^ 993 ' 1 8 i 80 7' 

27.200 0.240 1.328 0.966 0.580 7. '346 19.854 

" "■“2?,6U0 0.245' 1.302 0,958 0.591 ”o.6o5 20.935 

28.000 0.251 1.2/4 0,950 0.604 5.941 22,059 

’ 28.400 0.25 7"i.246 0.943 0. 6 20 "" 5 . io 7 23.233 

28.800 0.2o3 1.215 0.935 0.639 4.3J0 24-470 

" 29,200 0.2o'9 1,183 0.929 '0.6 63 " 3.413 2 5- 73/ 

29.600 O-2/o 1.14/ 0.925 0,696 2.3V4 27..206 

30.401 0.294 1-066 0.948 0.846 0.000 3“'),-f01 


Ml = 

= 2,600 


GAMMA 

, = 1.400 


SIGMA 2 

SIGMAS 

SIGMA4 

P2/P1 

P3/P2 

P4/P1 

<p4/pn4 


39,351 

57.863 

90.000 

3.004 

2.570 

7,720 

'4,850 • 


39,659 

57.091 

89.792 

3.04 6 

2.535 

7.720 

^.350 


40. 14 0 

57.949 

89.460 

3.111 

2 .-f-Oi 

7.719 

■ '4.851 


^jO.627 

58.029 

89,143 

3.177 

2.429 

7,710 

4, 851 


■4i;i2r 

'58.129 

'80.814" 

■~3;'24£“ 

"2 .'3 7 9" 

■ l~:j 17" 

, 8 3 


41.621 

58.251 

8 8.-^r62 

3-313 

2.329 

7.714 



42.129 

58.396 

■'88.14 6' 

"•3.382 

'2 .230 

7‘,7L2' 

4 .85 6" 


42.645 

58,565 

87.805 

3. 453 

2.232 

7.708 

4 .8 58 


43, 16 9 

58-760 

87.459 

“ 3,525 

2.106 

“7.704' 

■“4.860 


43. 701 

50.981 

87.107 

3.5 98 

2.U0 

7.700 

4,84=-'- 


'44'. 242 

59.231’ 

■ 

'“3,6/2' 

■2;095’ 

7:-695" 

^ . 8’.-:7 


.4^, 793 

59.511 

06,381 

3 . 74 8 

2.051 

7.609 

4.871 


45 ; 353 

59,824 

36.005 

■"3,825 

2,008 

■ /.6 02" 

4 .875 


45.924 

60.172 

85,619 

3,904 

1 .966 

7.674 

.881 


46. 507 

60.557 

05.221' 

■ 3,904 

1.924 

7.665 

4 » 0G6 


47.102 

60.985 

84-. 810 

4.066 

1.803 

7.655 

4.393 


'47. 7ir 

~" 6 1 : 457" 

’■84'. 3 85 ■ 

4'; 14 9” 

’1': 042" 

7':6'4 4"' 

" 4.9 90 


48.334 

61,981 

03,942 

4.235 

1.802 

7.. 632 

4,900 


4 8. 9"/ 3 

62,560 

83.479 

“■ 4.322 

1 .763 

'■"7.618 

“■"4.910 


49.629 

63,203 

02.993 

4.411 

1.7 24 

7 « 0 03 

4.928 


'50.305 

63.918 

B2.479' 

“4,503 

'1 .605 

""7.585- 

9:0" ■■ 


51,003 

t^r.717 

81.933 

4.597 

1.646 

7.565 

-'‘..95 3 


51.720 

65.6 13" 

"01.34 9" 

4'.6’94' 

T.'Bor 

7 .* 5 4 2 " 

4 .'96 9 " 


52.470 

66.626 

80,716 

4.794 

1 ,567 

7.5 15 

4,987 


53,263 

6 7. 781 

00.025" 

■4.898 

1.528 

"’*7.403" 

“■=.009’ 


54. 008 

69.1 13 

79.257 

5.007 

1 .4 0 7 

7,44 6 

__ 5 .03-4 


54. 962 

■ 70.675 

"70. 389" 

"■ "5, 12 0 

1.445 

“'"7.401 

5.055 


5 5 - 89 8 

72,547 

77.382 

5,?M 

1.401 

7,344 

5 ,105 


56. 914 

?4. 865 

76,169 

"5:370 

\ '1 1“ i- 

J, 4 ^ , 

■■■■/. 2 6 9 ' 

■9 , 1 5 r 


5 B . 0-; 5 

7 /.B06 

7s .62 2 

5.611 

1 .300 

7, L65 

5 .231 


o 1 . 004 

90.000 

68.422 

5 . 0c7 

1.134 

6.653 

■ 5 19' ■■ 

f ' 





TRIPLE POINT SOLUTION 


DELTA’ 

P2/P01 M2 . 

tAU 

.or.LTA3 0ELTA4 

SIGMA2 

SIGNA3 

S IGMA4 

P2/P1 

P3/P2 



P4/P1 

( P4/P11 


10.027 

0.123 

1.910 1.362 

6.438 

19.027 0.0 

38.329 

54.361 

90.000 

3.351 

2. 680 

8.980 

5.584 


10. <^00 

0.126 

1.891 1.160 


18.676 0.724 

33.770 

54.886 

89.728 

3.420 

2.626 

8.980 

3.584 


10. POO 

0o1 2’-5 

1.871 1,157 

0.488 

18.297 1.503 

39.249 

54.930 

89,435 

3.495 

2.5 69 

8,979 

5.5 84 


20,200 

0.1 32 

1.551 1.154 

0.4 89 

IT.q’5 2.285 

39.733 

54.903 

89.141 

3.571 

2.514 

8. 978 

5.535 


20.600 

0.1 34 

i.«.'’0 1.3 50 

0.490 

17.529 3,071 

40.223 

55.075 

.'8 3.844" 

3 .648 

2.461 

8,976 

5.586 


2; .000 

0.1.2 7 

1.3,10 1.‘'46 

0.490 

] 7. 140„3.8 60 

.40.71 9 

55.177 

88.545 

3.726 

-2.409 

3.974 

3.58 3 


21.400 

0.140 

1.789. 1,142 

6.401 

16.74'7 4.6 5 3 

41.223 

55, "^00 

88.243 

3.305 

2.358 

8,971 

5. 590 


21 , 8 00 

- 0 . 1 4.3 

1 .768‘. 1.137 

.0,493 

1 5. 350 5.450 

.41.733 

55.444 

-87.937.. 

3 .386 

-2.308 

8.963- 

5.592 


22-200 

0.0 46 

l.f67' 1.332 

6.404 

15. '94 8 6.252 

42.250 

55.611 

87.628 

3.968 

2.259 

8.964 

5.595 


22.600 

0,1 ‘*'1 

1.726 1.127 

0,406 

15,541 7.059 

42.775 

55.30’ 

87.314 

4.052 

? . ? 1 1 

fl-QAO A AQO 


2? .000 

0.1152 

1 . 705 • 1. ■< 21 

0.497 

15.12 3 7,8 72 

4J.307 

56,01 f 

86.994 

4.137 

.-Vi. ^ > ., -8„ 

2.165 

8.955 

. 5,60? 


2 3 • *T 0 0 

0.156 

1.683 1.135 

0.409 

1,4. 71 3.6 90 

43.843 

.56.259 

v86.66Q 

-4.223 

2.119 

8.949 

3,606 


23.800 

0.159 

1.662 1,109 

0.502 

14.284 9,516 

44.399 

56.529 

86.337 

• 4.3U 

2.0 75 

8.943 

5.611 


2i.200 

0.162 

1.640 : .303 

_0. 504.. 

13.851 1 0 . 3 -^'9 

. 44,957 

56.8?8 

85. 998 

4.400 

2 . 031 . 

a. 935 , 

5.616 


2 - , 0 0 0 

0.165 

1.618 1,096 

6.507 

13.410 U.190 

45,526 

57.160 

85.650 

4.491 

1 . <^88 

8.927 

5.622 


23.000 

0. ■’6*’ 

■> . 595 1 .088 

O.'^IO 

12.960 12.040 

46.105 

57.525 

85,292 

4.5 83 

1 . ^46 

8 ,9’ 8 

5 . 6’ 3 


28.^00 

0.172 

1 -572 1 .081 

0.514 

12.500 12, .900 

46,695 

57.928 

34.923 

4.077 

1.905 

8.903 

5,636 


- 23. POO. 

,0. 176„ 

1.549 1.073 

-.0.518. 

,12.030. 1 3,770. 

,47.298. 

58.372- 

84.542_ 

-4.773- 

-1.864, 

E.897-- 

5.64 4_ 


2 h . 2 0 0 

0.179 

1.526 1.065 

0.52 2 

11.547 14.653 

47.914 

58.860 

84.147 

A .871 

1.8 24 

8,885 

5.653 


26.600 

0.1P3._J .503_1.05 7_ 

-0,527.. 

i:.051 15.540 

. 43.545 

39.397, 

33.736- 

-4.971- 

-1.7-8 5 

8.3 71- 

5,663 


27.000 

0.187 

3 .478 1 . 048 

0.532 

10.539 16.461 

49.191 

50.990 

83.306 

5.073 

1.746 

8.886 

5.674 


-J!Z_*4.00. 

0 ; 1 ’ 

1.6^4 ’.019 

0.538 

iQ..0l.U.:j,7.3 8o 

49 * 8.54 

60. 644 

8’. 8 55 5-178 1 7f)7 

A AQ7 


27.P00 

0.3 95 

1.429 1.030 

0.545 

9.463 18.337 

50.5,37 

61.369 

82.378 

5.285 

1,669 

8.819 

5,701 


2P.?00_ 

0.1 SP. 

-1..4.03_i.»020. 

-0.5 5Z_ 

_._8. 8 93 , 19 . 307 ,. 

51.242 

62.175 

.81.872 5,395.JL,633 

_3.797._ 

—5.718 


28.600 

0 ,203 

1.377 1.010 

0.561 

8.298 20.302 

51.971 

63.077 

81 .329 

5.509 

1.592 

■ 8.7 72 

5.7’6 


. 2«.000 

.0.207, 

-1.3,50.2.000. 

0.570- 

.-,7.673 ,21.327 

-52.729- 

64.091 . 

80.742- 

-,5 .626- 

-1.5 54 

-3.743- 

.5.758 


2‘5.6-00 

0.212 

1,3 2.2 0.090 

0.501 

7.014 22.386 

53.521 

65.P44 

30.101 

5.747 

1.516 

8.710 

5,783 


8 0.0 _ 

0.216 

JL,.2.a3 

0,5 04 

6.3’ 3 73.437 . 

54. .352 

66,5 68 

79-3 91 5-873 1.8,75 Q.ATH 

5 P‘ 3 


30.2 00 

0.221 

1.263 0,969 

0.609 

5.561 24.639 

55.232 

68.113 

78.589 

6.006 

1.436 

8.622 

5.35 0 


?0.600 . 

0.226_ 

,1 .221_0.959_ 

,Q.627_ 

_4.,746 _25.854,. 

56.175, 

69.95 ? 

77.662— 

-6.146- 

-1.39 3—- 

-8.562- 

-5,896 


33.000 

0 . 2 ■ 2 

1.196 0,049 

0.6 5 0 

3.850 27.150 

57.199 

72.227 

76.549 

6.296 

1.348 

8.48 5 

5.956 


2 1,400 

0.238 

1 ,159. . 0.041 

0.6P2-_ 

_2.845 23.555- 

53.33-3 

75,165 

,7 5.1 4 0,_ 

-6.460- 

1.297 

-8.373 . 

-6.0^1- 


3 i o 00 

0.26-5 

I . 1 3. 6 0.93 d 

0.7 28 

1.689 30,111 

59.660 

79.327 

73.166 

6 • 64 6 

1.236 

3.213 

6.175 



X,.O4,5__0,Q4ti._ 

„Q.,264„ 

_-0.QaQ-_32.3Q3_ 

61.8 89,. 

90,000 

6 7,.,99Q„ 

16.949- 

-1.107-:i_ 

— 7 ,6 95- 

— 6.619— p 

K) 


3.0 


Tal)le_9 ___Triple-Poirit_ Solutions 



TRIPl 8 

POINT 

sni.UT ir-M 


Ml 

= 3.000 


gamma 

= 1-400 


PELTA7 

P2/P0? 


M3 

M4 

0'3LTA3 

DBLTA4 

SIGMA2 

SIG0A.3 

SIGMA4 

_-P?/P.L- 

P3/P2. 

..P4/.PI,- 

„<P.4/P1)^1. 

■ l?.o5A 

b.’lO: 

2.01 7 

,1.743 

0.675 

1. 9.656 

0.0 

37.359 

5.2,319 

90.000 

3.700 

2.793 

10,333 

6,372 

20 » ago 

0.106 

1,973 

1 . 7 ? 6 

0.475 

1 3.97‘6 

1 .424 

33.238 

52.356 

89,501 

5,856 

2.6 80 

10,.3?3 

,..6'. 37 2 

20. BOO 

0.107 

1.951 

1 .2 .3 7. 

0 . 4 76 

.1 8. 60S 

2,192 

38,718 

52,400 

89. 232 

3.941 

2,621 

10,331 

6.373 

2' .200 

0. ) 10 

?. .0 20 

1.228 

0.4 76 

1 8.2 3S 

2 .9 6 2 

30,204 

52.463 

88.«61 

4,024 

„2..?64 

..J .0.130,. 

_6, 374 

21.600 

O.V 2 

1 .908 

1 .222 

0.4 77 

17. 866 

3.734 

39.695 

52.543 

88.648 

4', 117 

2.5 09 

10.328 

6.376 

22.000 

0 . U 5 

1.886 

1 .718 

0.478 , 

17.490 

-4.510 

40.192 

52.642 

.86,414. 

..4.206 

2i455 

10.325 

„ 6.378 

?2;AO0 

0.117 

1 .864 

1. . 2 1 3 

0.479 

1 7.1 1 2 

5.238 

40. 695 

52.761 

8,4.136 

4.2-97 

2.A02 

10,322 

6.380 

22.800 

0.120 

I .84 2 


o,.4yo 

16.731 

„ 6,069 

43 . 205 

■"2.900 

87.456 

..4, ??0. 

2.351 

.10,319. 

..6,383 „ 

? ’3.2 0(1 

0 . f 2 2 

1 .81 9 

1.201 

0.482 

16.346 

6,854 

41.721 

53.060 

87,572 

4,434 

2.300 

10,314 

6.3S6 

23-600 

0.125 

1 .707 

1.10 4. 

0.4p7 

1 5.'957 

7.643 

4?, 245 

53.242_ 

S7. 245 

4.579 

.2,251 . 

- 10-310.. 

__.6,?90 

2A.OOO 

0.127 

1 .774 

1 . i H 8 

0.486 

: 6.66? 

3.437 

42.775 

53 .445 

36 . 993 

4.676 

2.204 

10.304 

6,3 94 

?A.^O0 

0.170 

1,752 

1.180 

0.4 87 

15.164 

„ 9.2 3 6 

43.313 

53.675 

8 6.696, 

_.4.774, 

2.157 

1 0.298 

., 6.399 

?'-.B0b 

0.133 

3 ,7.?o 

1. . • 7 3 

0.489 

1 4.760 

10.040 

A3. 860 

53.930 

86.393 

4,874 

2 , .1 1 1 

10,292 

6.4 04 

23.200 

0 . 7 ^ 

] .706 

1 .165 

0 . 49 7 

14.350 

1 0,35 0 

4,. ^ 4* I 

54.211 

86 i 044 

-4,976 

2,06'7 

10. 2?4.. 

.- 6. AlO 

25.600 

0. 13>7 

1 .68.? 

1.157 

0.495 

13.933 

11 .667 

44. 973 

54.521 

35.768 

5.079 

2.023 

10.276 

6.41 7 

26.000 

0.1 

1.65 9 

1.143 

0.493 

>3.508 

1 2.4 92 

'45,551 

54,862 

35.444 

■",144 

J. -98,0_ 

..10. 2.6 7.. 

6,424 

2 6.^ 0 0 

0. 144 

1 .635 

1.140 

0,5 0. 

13.075 

13.325 

46,135 

55.236 

85.1 n 

5 . 29 1 

1 . 939 

10.257 

6,432 

? D . s 00 

0.147 

1 .61 1 

1.13 0 

0.6 06 

12. 634 

14.166, 

45.729 

55.646 

,84. 75 g 

... 5,4 00 

1.397 

JO. 2 46. 

6,440 

27.200 

0.150 

1 .587 

} . 1, 2 1 

0.609 

12.182 

15.0,18 

47.335 

56.094 

84,414 

5.5> 1 

1.857 

10.234 

6.-450 

77.600 

0.167 

1.567. 

?. .,1 1 1 . 

.0 , 5 J. 3. 

n.719 

15.831 

47.954 

56,586 

84,04 7. 

. .5.6 24 

1,81 7. 

10.220- 

6,461 „.... 

73.000 

0.156 

1 .537 

1 . 1 0 i 

0,618 

1 1 . 2‘4 3 

1 6.757 

43.586 

57.1 24 

8 3.6 64 

5,739 

1 ,778 

> 0. ?05 

6,47.3 

7^ .AGO 

0.‘i 

1 .51? 

1^00 i. 

0-5 2 3. 

10, 753 

.1.7, 6 47.. 

49.2’4 

57,715. 

.83,. 26 5.. 

5.3-5 6.. 

.1,740.., 

...1 0,.1.33. 

6 .436 

2 3, BOO 

0.163 

1 .486 

1.080 

0,6 20 

10.747 

■ 8.563 

4«. 399 

58.165 

82 * 846 

5.977 

1.702 

10. 1 70 

6,501 

27.200 

0.;a6 

1.460. 

1.06 9. 

0-535 

...9., 72 3 

.] 9.477 

.,50. 5 52 

59 . 0 82 

82.404 

..6.100 

1.664 

10.150 

.. 6.51 7 

27.600 

0.160 

i .47^ 

• .053 

0-543 

9.179 

20.421 

51 .2 36 

59.876 

S1-.935 

6.226 

1.6 26 

lO.'l 27 

6,536 

3 0 , \)00 

0,;73_ 

i .406 

.J. .04 6. 

0.5 51 

. 8,610 

21 ..390 

52 . 014 

.60,7.59. 

.31,434 

„6.256 

.1 .587 

10.100. 

._-6,557 

30.A00 

0.177 

1.378 

I .07 4 

0.660 

8.014 

2 2., 3 86 

52.7o9 

61.747 

B0..49*t 

6,490 

1.552 

10. '3 70- 

6.5 82 

70. 700 

0.1 HO 

1 , 7 Q 

I.o2.2_ 

0..5.9 0. 

...Z. 38 6 

.23.414.. 

. 6 ? . 5 5 0 

.62.864_ 

.30.306. 

__.6.6 2B. 

1-514 

..2 0,0 36. 

_6.610 

3i . 700 

0.18.. 

t .31 0 

: .01 •) 

0 , 5 8 2 

0.71 9 

24.^41 

54,381 

64.139 

79.654 

6 i 77.2 

1.47 6 

. 9,995 

6.644 

.600 

0 . : 8 3 

1 . 7. 8 7 

0.097 

O.50? 

6,004 

2 5.5 96 

35,253 

6 5.61 5 

78.930 

..6.72.2. 

.1,437. 

9.946 

6.6.84 

7 2.000 


1-26 6- 

0. 18 4 

G - ( 1 4 

5.700 

2 6 . 7 7 0 

56.1 .8.?. 

6 7. '15 8 

78.096 

7.091 

1,39 6 

9.887 

6-734 

3 7.^00 

0 

• •» o 

i * *- , . 

0.07 1 

0 .635 

4,320 

2 ? . 0 3 0 

57.1 87 

to .474 

77,109 

. 7., 2 SO 

1 . 3 ■" 3 

9 . 8 1 1 

. fc , 799 

' 1 . 8 0 ) 

0 ..-07 

1 . ' 7 .) 

0. >•' 7 

0 . /■ 6 3 

3 , 4 2 .8 

2 >.'''■7? 

'= .3 . P 9 ■> 

72.159 

7A.r44 

7.4 34 

! . 3 0 6 

9 . 7 09 

6,3 8 6 

7 ri'o 

0-. 7 08 

^ ^ *1 *7 

*, t , . t 

0 . 0 4 •} 

0.7 02 

2 - 3 ^ 

3 0.-370 

5 9 . 5 7 .1 

7 5 , f 2 5 

74-235 

. 7.640 

1 ^ ’ 

9.5 b 8 

.7.014 . 

"At , 

* • 

0 . ? : ■ 

1 .0 7? 

.0 . 0 h 6 

0 - .7 8 0 

0.000 

3 ( . 8 7 > 

62.607 

90.000 

67.712 

ci *11 

1 , 

r .823 

7.713 H 



Table 10 Triple -Point Solutions, = 3*-2 


TRIPLF POINT SnulTION ' ' Ml = 3.?00 GAMMA = 1.400 



nn.TA 2 

P2/P01 


M3 

M4 

DFLTA3 DSLTA4 

_Sir,MA 2 

Sir,MA3 

SIGMA4 

.-P 2 /P 1 

P3/P2 

T V V 

* 

, P4/.PJ 

_LP4/P1J'* 


20,1 00 

0 . 0 B 2 

2.115 

3.32 1 

0.464 

20.100 0.0 

36.451 

50.145 

90.000 

4.051 

2,908 

11.780 

7,214 


20.500 

0. 034 

2 . 0^2 

1.317 

0.4 64 

19.747 0.753 

36.919 

50,146 

-.5^*751, _ 

,4.3 44., 

.2.843 

—Il.T.OO. 

—7,214 


21.000 

0.036 

2 .064 

1.312 

0.465 

19.304 1.696 

37.511 

.50.171 

89.440 

4.26 3 

2.763 

11.779 

7,215 


21.500 

0.03P 

2.0^5 

1.306 

0.465 

13.35 9 .2.641 

38. Lil 

50.220 

89.127 

4.3S4 

2.686 

_Jl-l *,7.7 7_ 

7*216 


2 2.000 

0.091 

2 .006 

1.29 9 

0.466 

18.432 3.588 

38.719 

50.297 

88.81.3 

4.507 

2.612 

11.775 

7.213 


2? . 500 

0-004 

1.977 

1.292 

0.467 

17.963 4.537 

39.335 

50.4 00 

..88.496 

.4,6 3.3. 

,2.541 

-1 1,772. 

7.221 


25.000 

0 . 0^6 

i .94)5 

! , 234 

0.468 

17,510 5.490 

39.960 

50.531 

83.176 

4.761 

2.472 

U. 7 68 

7.224 


25.500 

O.O'^O 

1.919 

1.276 

0.470 

17.054 ...,6.446 

40. 594 

50.692 

,e7.s852._ 

-,4,892. 

..2,40 5 

__11 .763, 

-7,228 


24.000 

0 - 10 ? 

I . 3 .3 9 

1,267 

6.471 

16.593 7.407 

41.238 

50.883 

87,524 

5,024 

2.3^0 

11.758 

7,233 


24.500 

0.1 04 

1,360 

3. . 2 5 7 

0.474 

16,127 8.373 

4). .892 

5 1 . .1 06 

87*19 0 5.. 1.60. 

-2.27,7' 11,75 L- 

7.238 


25.000 

O.IO^ 

1.350 

1.247 

0.476 

15,654 9.346 

’42. 556 

51 .363 

B5,85). 

5.298 

2.217 

11.744 

7.244 


2 5.5 00. 

.o-"! 10 

.1 *300_ 

-1.-.Z3 6 

0.478 

1 5. .175 .1.0.325. 

.i^3.232. 

51.656 

8,6 , 504_ 

-5.438 

2.15B 

__11 .736. 

„7*2 52.._ 


25.000 

0.113 

1 .769 

1.22 5 

0.4 Ri 

14.589 11.311 

43. 920 

51 .987 

86,149 

5.582 

2.101 

11.726 

7.260 


25,500 

0.316 

1 .7 39 

i . 2 1 3 

0.48 5 

14.193 12.307 

.44*621 

..52.3 59 

..85*784 

-,5,728_ 

.2,045 

11.715 

,_7.269 


27.000 

0.119 

1 .70.3 

1.201 

0.489 

13,687 13.313 

45.335 

52.776 

85,409 

5.877 

1.992 

11 .703 

7.279 


27.500 

0 . 1 3 

1.67 7 

1.13 8 

0i498 

13.170 14.330 

46.065 

53.241 

e5,,a.2i 

6.0,39 1,939 

- U .590 

7.290 


2B.000 

0.125 

1 . 6 4 

1.175 

0.497 

12.640 1 5.360 

46,81. 1 

53.760 

84,618 

6,184 

1.888 

11 .6 75 

7.303 


2"».500 

0. 1 2=> 

1 . .6 1 3.. 

1.161. 

.0,50.3.. 

.12. ,095 .16,405 

.4-7.574 

54.337 

.84.19 8..„ 

. 6 . 343 . 

a. 830. 

1 1 .65 8-. 

—7.317 


2^.000 

0.-152 

.1.581 

i'. : 4 7 

0.5 OS 

11.532 17.468 

48.353 

54.981 

83.759 

6.505 

1.789 

11.639 

7.334 


20.500 

0 . 1 3 5 

i.,5.4 3_ 

.1 . 1 32 0.515 

10..950..18.550 

49.163. 

,5 5.7 00 

,83.296. 

,6.672. 

-.1.741. 

11.617,- 

.7.35? 


30.000 

0 . 3. 3 a 

1.514 

1.317 

0.5 22 

1 0.344 19.656 

49.994 

56,50.6 

82.806 

6 .84.3 

1.694 

11 . 593 

7.373 


30.500 

0.14? 

1 .4 80 

1. U)1 

0,530 

9.71 2 2 0^788 

.50.8 52. 

,,5.7.412. 

,J2..2 82™ 

„X.OIS_ 

.1,,.64 8 

--11..5 65 

.7,.39a 


31.000 

0 . } 46 

i .'+45 

1.085 

0.539 

9.043 21.952 

51.744 

58.439 

31.717 

7.200 

1.602 

1 1,532 

7.426 



31.500 

0,3.40 

.1 .-0? 

l..i36 3. 

.0.5 50 

,,8.346_23.154 

52,674. 

.59.611 

,.BX . 1 02 

-.7.283, 

-1,556. 

—11.494-. 

—7.459 


32.000 

0 . : 5 3 

1 .371 

1.050 

,0.563 

7.597 24.403 

53.651 

60.964 

80,421 

7.583 

1.510 

11,449 

7.499 


32.500 

.0.15 9. 

1 . ■' 3 2 1 . 0 33 

0.579 

„i>. 79 0.25.7 10 

.54. 633. 

.62.5 52 

..79.6 54- 

7.788.. 

-1,46 3- 

11.395- 

7,54 7.. 


33.000 

0.162 

i.391 

1.014 

0.5 96 

5.910 27.090 

55.793 

64.457 

78.768 

8.004 

1 .<^15 

11.327 

7.608 


33.500 

1 4 7.. 

.1 . 24 0 0,996 .0 .63 9 

_ 4. 92.9 2S.571. 

56.999 

,6 6.319 

-77.704 

.,3.23 6., 

-1,3 65- 

n.,.233 

7.6 38 


?4.000 

0.172 

1 . 1 9 s 

0.Q77 

0.650 

3.807 30,193 

5B.350 

69.909 

76.352 

8.491 

1 .309 

11.115 

7.801 

— < 

34.500 

0.173 

1.142, 

.0.960 

0.693. 

„2.46 4 3 2,0‘3 6 

59. 940 

76.3 7 7 

,74.433 — 

-.8.783 

-1.243. 

.-10.920 - 

—7.936 


35.1 HO 

0.131 

1 .0^1 

0.970 

0.593 

O.OnO 35.180 

63.199 

90.000 

67.534 

9.351 

,1.07? 

10.0 36 

• 8.897 


h-" 

+T?* 

4 ? 





Table 11 Tri ple -Point Solutions, M-, = M'.O 


TPIPLn POINT SniUTIQ.NI Ml = 4.000 GAMMA =M«4no 


nSLTAZ 

P2/P01 


M2 




0PLTA3 

D8LTA4 

SIGMA? 

SIGMA3. 

..S.I.GMA4. 

P2/P.1 

.P.3/P2. 

_.__P4/..P1.__ 

..I.P.4/.PDTL. 

2 0,8'i4 

0.036 

2 

,507 

1 

.6X0 

0.455 

20.854 

0.0 

33 .444 

43.921 

90.000 

5.503 

3.362 

13.5 bo 

11.130 

21.600 

0.0 33 


.4 53 

1 

. 598 

0.435 

20.265 

1.335 

•^4.314 

M.892 

.8 9.526 

5.765 

_3,.2p9. 

—18.499. 

._U.131__ 

22.200 

0.030 

2 

.41.0 

1 

. 587 

0.436 

19. 794' 

2.406 

35.023 

43.904 

89.326 

5.982 

3.092 

18,497 

11,132 

22. POO 

0.041 

2 

.3 67 

1 

.574 

0.4 36 

19.3 2-5 

3,475 

35.742 

43. 947 

89i 025 

6 . 203 

2*_?-82- 

_1.1- 4-9 5_ 

.1.1.. 1 35 

23.^00 

0,0-'- 2 

2 

. ^24 

V, 

.561 

0.437 

18,856 

4.544 

36.470 

44,023 

.88.722 

6.42 8 

2.876 

18,^91 

11.139 

24.000 

0.04 6 

7 

.28 J 

1 

. 547 

0.430 

1 

5.612 

3 7. 207 

44.131, 

8 8.41.8 

6.659, 

2.776,. 

__18.486. 

.U-144__. 

2 '^.6 00 

0 , 0 '* 3 

c. 

.2 35 

1 

.5 32 

0.44 0 

17.91 8 

6.582 

37.955 

44.272 

.38.111 

6. 8«4 

2.680 

18.480 

11,150 

2=^. 2 00 

0.0‘-Vl^ 

2 

. I o? 

X 

.5X6 

0.442 

.17,447. 

.7.753. 

.3 8.713, 

44,448. 

.87.301. 

_7.135. 

_.2 . 5 8 9.. 

_„I8,4 7i. 

.11,159 

2?. 800 

0.040 

2 

. 1 5 2 

I 

.500 

0.4A4 

16.973 

8.327 

39.481 

44.659 

87.487 

7.380 

2.502 

18.464 

11.166 

26. -^00 

0.030 

2 

. loo 

3 

.482 

0.446 

3 6. 496 

9.904 

40.2ol 

44,908 

87 i 168 


._2..41.9„ 

— l.a.».45.4_ 

.iL.lI6 

27.000 

0 .052 

2 

.046 

X 

.464 

0.449 

1 6 . 0 L 4 

10.986 

41 ,053 

45,196 

86,843 

7.385 

2,339 

18 .-4^-3 

U . 1 8 7 

2 7.600 

0 .054 

? 

.02? 

1 

.445 

0.452 

15, 526 

12.074 

41.858 

.•43,524. 

e6,510_ 

_.a-145. 

..2,263 

—.18, 431 

.11-200, 

23.2 00 

0.05 5 

X 

.7 79 

1 

. 425 

0.455 

15.031 

13.169 

42.676 

45,897 

86.170, 

8,410 

2,190 

18 .417 

11.214 

?S. f^OU 

0 .05,7 

1 

, « 3 c 

1 

.404 

0.459 

1 4,5?B_ 

1.4.272 

43.509. 

46.317 

_85.820 

__8«681 

2.1 20. 

_..l 8 ,401. 

.11-230 

20. 

0.059 

i 

.50 2 

1 

.383 

0.463 

3 4. 015 

15.385 

44.358 

46.737 

85.459 

8,958 

2.052 

18.383 

11.249 

30,000 

0,06’ 

1 

,040 

1 

.361 

0.46 8 

13.490 

16.510 

45.224 

47.314 

85,085 

_9.24 0_ 

..A=.,a87_ 

__J 8.3.61- 

7 1.. 2 69 

30,600 

0,063 

1 

4 . 

.9 05 

1 

.337 

0.473 

12.951 

37.649 

46, 109 

47,901 

84,696 

9.528 

1.925 

18.341 

11.292 

3 : . 2 0 0 

0 . 065 

X 

.7o0 

1 

.316 

0 . 479 .. 

.12.397. 

1 8.803 

..47,0X4 

48,5 57.. 

34 , 290. 

_.9.32 2, 

.JL,B65_ 

18.315 

.11.318 

3.1 .800 

0.067 

1 

.73 6 

1 

.292 

0.485 

11.823 

3 9.977 

47.943 

49,290 

83.863 

10.124 

1,806 

18.287 

11.34B 

32.^00 

0-060 

J 

. 6 7 0_ 

T 

_ t 

. 7.6 7 

.-Q.4,9 2_ 

.,11.228. 

21 .172 

..48,693. 

30.111 

_R3,4U. 

.10.433, 

.1,750. 

— la.254_ 

.11.36 1 

33.000 

0.071 

1 

.62 5 

*> 

X 

..•’42 

0.503 

10,607 

2 2.3 93 

49.B83 

51.034 

82.930 

10.750 

1.695 

3 S.?17 

11.420 

3 3, ^-00 

0.073 

jl 

.578 

1 

.216 

0. 510 

9.955 28.645 

_.50..903.. 

.52,07.8. 

_82.4i.2_ 


.1.641 

_.1.8 .J..T5_ 

.11.464 

^4.200 

0.075 

1 

.5 30 

1 

.15 7 

0.520 

9.265 

24.935 

51.964 

53.269 

81.850 

11.413 

I.5BB 

IR.l 25 

11.51 7 

3-^. POO 

0.077. 


,.4 81., 

.U162 

„0.?33_ 

-_a. 5 23.. 

,26.272 

53.074 

54-641 

..81 .230..11.762- 

..1,536. 

._18-066. 

.11,579 

35.^00 

0.050 

? 

,431 

.1 

.13 3 

0 . 547 ' 

7.73 3 

27.667 

54.245 

56.247 

80.534 

12.1 27 

1,434 

17.995 

11.655 

36.000, 

0,052- 

1 

,3.71. 

1 

.103. 

_0 » 565- 

__6.3 6 3.. 

29.137 

. 55 . 495 .. 

.58,166 

. 79.735 

12,510. 

_1,4 31 

— 17.907- 

-11,750 

3 6 , 6 0 0 

0.055 

1 

.5 2 3 

1 

.073 

0.58 7 

5.839 

30.71 1 

•56.852 

60,533 

78.784 

12.919 

1.377 

17.794 

11,874 

3 7.20 0 

0.05 5 

1 

•J 6 0- 

•» 


Q.. 63 J 

4. -7 .4.6 

J2..43 4„ 

38. 36 5.. 

.63,6 06 

77.537 

.13. .3 6.5. 

_1.320 1.7.6 38- 

_12,04 7. 

37.R00 

0 . 09 1 

1 

.loo 

X 

.00 6 

0.661 

3.397 

34.403 

60.138 

67.989 

75.921 

13.872 

1,254 

17 ,395 

I 2,324 

3 3.722. 

0.100_ 

-X 

.01.1. 

.0 

, .93i. 

.0..,929_ 

,._0.00 0 . 

3 8.722. 

64.77.7, 

90.000 

.67.295 

..15,110 

-1.040. 

._1 5 . 7 3 9 . 

.14.497— 
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Table 12_ Triple-Point Solutions, = 




T6 IP1.P 

PHI NT 

SOLUTION 


Ml 

= 5,000 


GAMMA = 1.400 




P2/P0 

1 ^2 

•• M8 

M4 

delta? 

DELTA4 

.SIGMA 2 

S lOMA? 

SIGMA4 

.. P2/P1 

P3/P2.. 

_P4/P1_. 

{ P4/.PU4 


2 0.862 

0 . 0 1 -i- 

2.64^ 

t.=io 

0.415 

20.362 

0.0 

30.796 

39.271 

90.000 

7.479. 

3.87 8 

29,000 

1 7,2 50 


?’ .600 

0 .01 5 

? .8 7 3 

' .897 

0.4] 5 

,70.344 

L..256 

31.659 

3'2. 733 

89,686 

7.86 8 

3.6R6-. 

28.9G9 

17,251._,. 


2 ? . 2 00 

0.015 

2 .818 

1.377 

O.'tlS 

1-9.02 3 

,2.272 

32.368 

39.234 

89,431 

8.103 

3,539 

28.997 

17,254 


• 72.600 

u.o: 6 

2.763 

1 . Ht.0 

0. 41 6 

to- 5-1 6- 

3.734 

33.085 

3 ’,2 63 

89,1 77 

a. 57c; 

3.401 

„?8 .99.4. 

17.257 


20,^00 

0.017 

2.70 3 

1 . 8 4 2 

0.41 7 

19.107 

4.293 

.^3. BiO 

39.520 

8 3. 9>..’ 

8.864 

3.270 

28.990 

17.262 


2'-. 000 

0.01 7 

2.658 

1.823 

0.41 8 

1 8. 700 

5.300 

34,542 

39 .405 

88.668 

9,210 

.’,1.47 

28.984 

17,2 68 


?-'‘.6O0 

0 . 0 J 2 

2.50' 

1.802 

0,414 

I B. 704 

6.306 

3 5.. 7 82 

39.518 

8 8. 4U 

9.564 

3.030 ' 

28,978 

17.276 


25.200 

0,0* ’ 

2.5-8 

T .78 1 

0.491 

17- 389 

- .7.311 

36.030 

?'7.66 0 

88,153 

9.9 2 5 

2.919 

.20.970 

17.28 5 


2‘=5.800 

0.016 

2.49 5 

' . 758 

0.422 

17. 484 

3.31 6 

36.786 

’ 9 i r .3 0 

87.8°2 

10.2’2 

2o814 

23.961 

17.295 


?6,^00 

0.020 

? .44'3 

1.737 

0.4 74 

] 7.077 

9-372 

’7.550 

40i0’0 

87.6-2 9 

} 0,657 

2. 7"' 4 

_2-8..9 50^ 

.Z.7.3 01L. 


" 7.000 

n.o’i 

2.351 

1 .7X 3 

0,426 

16.669 

10.331 

?S. 324 

40 i 261 

67.362 

1.1 .04’ 

2.63 9 

28,938 

17.’2i 


27.400 

0.022 

2 .3’"> 

1 .68 9 

0.42Q 

1 6 . .7 5 9 

i 1,341 

’9. 107 

40^523 

87.071 

n .^’8 

?..5?9„. 

. 28.9,2 5 

1 7.336...: 


28.200 

0.022 

2 , ? H 8 

1.66 4 

0.431 

: 5. 345 

12.355 

39 . 900 

40.018 

B6.315 

11 . B 5 4 

2-A43 

23.910 

1 7.353 


2S.B00 

0.02’ 

2,737 

■ ,639 

_0.-34 

. Z:>* ^26 

1 3 . 3 74 

4-0.704 

41 ,147 

.86. 5 ’4 

.12.238 . 

2.3 61 .. 

.20.093 

17.373. 


2^.600 

0 .0?4 

7.1 56 

1.61? 

0.437 

!, 5.002 

14.393 

41. 51B 

41.514 

86.246 

12.648 

2. 28? 

2 8.075 

17,394 


’0-000 

0.02 5 

2.136 

1.8 36 

0.441 

1 4, 571 

15.420 

42.344 

4 1 , 9 !, 8 

85.950 

_Z3*.967_ 

_2.2 08 

20.. 8 5 5. 

.1 7.41 7 , . . 


20.600 

0.076 

2 . f) 8 5 

1.559 

0, ^ lif 4- 

14.1 3? 

16.4,57 

47.183 

42.365 

8 5 . 64 7 

13.402 

2.137 

28.832 

17 . 444 


31.200 

0 .0^-6 

2. . 0 7 5 

1.537. 

.0.449. 

_i 3.6 36 

.17.5 14 

44.036 

42.856. 

85. 333 

..3 3. .92 6. 

.2-069..„ 

_23.307 

17.473 


2 i - 8 00 

0-027 

1 . vB5 

1 .503 

0.453 

13-2?. 3 

1 3.57 7 

^ ^ ? 0 ^ 

43.396 

85.009 

14.368 

2.003 

2S.779 

1 7.505 


^2.^00. 

0.02 0 

J .'-'iS. 

_:.475„ 

..0,4 5 8 

1 2 . 759. 

19.641 

45.7BS 

43. ’39 

84.672 

14.812 

1.940 „ 

..23,743. 

17-541 . 


’7.000 

0-02? 

1 .8 85 

1 .44 6 

0.4 6^ 

12.275 

2 0.7,2-3 

46.600 

44.647 

84.320 

15.276 

1.8 80 

28.714 

17.581 


23.600 

0.030 

l.8’ = 

1 . 41 6 

G.470 

U»?.76. 

.,2 1,8 2 'Z 

47.6,12 

.45.3 61 

.d3...,9.5,l„ 

...I5..744 .. 

-1.3 21 

.2 3.676.. 

.17,626 


2.00 

0.031 

1.734 

1.336 

0.'-77 

11,258 

22.942 

48 . 556 

46.155 

83.562 

16.22? 

] .765 

28.633 

1 7,6^6 


24.000 

0.032 

1 ,;73 J ._ 

_1 .35 5- 

.0.4 35 

..10.719 

24.031 

49.526 

47.036. 

a.3 - 1 5 0 

16.711 

1 . 7 1 1... 

.2 3-. 5 85 

17-73 3.. ... 


.60 0 

0.033 

\ . 6 8 ? 

1.324 

0.'‘93 

10. 1.54 

25,246 

50.525 

48.018 

82.709 

17.7] 2 

■1 . 6 5 B 

28.530 

17.799 


26,000 

0.0 ’4. 

..Z .6,3 0. 

.1.291. 

0.50 a. 

. 9.558 

26.442 

.,51. 558 

49.119 

a2.2'34 

17.726.- 

1.606... 

..2.8,467. 

17.874 


2 6 . 6. 0 C 

0.035 

I . 5 7 8 

1 .259 

0.5] 4 

8.0,7 6 

27.674 

52.632 

50.366 

81.716 

18.256 

1.555 

?8, 395 

1 7.961 


2,7. 20 0 

0.0 34 

1 .5 24 

i-...7 2.‘> 


8 . 7 40 

.28.951_ 

.53.754 

3.1.7 94 

.,8 1.14 5, 

-ia.804_ 

.1 .505... 

„2a.309. 

18..065._ 


3 7.80 0 

0.03 7 

1.468 

1.190 

0.542 

7.515 

30.2fl5 

54.938 

53.455 

00.50? 

19,37-5 

1.4 56 

20.206 

18.1 92 


..■2 3.-00 

Q.O? 3. 


.l.Z5 4„ 

-.0.5.60- 

„ 6.7 08 

.31,6 92. 

.56.201 

.55,430. 

-79.762. 

.1 9.97.4... 

1 .4 06- 

-23.079 

.IB.,349 


37.000 

0.039 

1.34a 

1 . U 7 

0-5 82 

5.300 

33,200 

57.571 

57.856 

78.879 

20.613 

1.354 

27.915 

18.534 


36.600 

_0.i]4Q. 

’.782 1 .077 

_0,6 12 4,744 

-34.856_ 

_59,iOI 

..60,9 97 

7.7..-.765 7L.30a_i 2.99. 2.7^690. 18-841 


40.200 

0.042 

1.206 

1.035 

0.6 57 

3.441 

36.759 

60.902 

65.634 

76.206 

22.102 

1.237 

27.342 

19.299 

— 


JL.-.0J..CL_O. ,.99.a__a,53A_.0 *-0i)a_6tl . Q.9.?„ 

..65., 823. 

_.90.000_ 

.6.1.32 0. 

.24.i0S_JL-0 2.’_. 

_,24.666_ 

-23,.3 35— K- 
0 
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Table 13 ^Triple-Poin-^ Soj^tio^s , 



• Tkl'PL*^ 

pni'NT- 

sn'L'j7. 

IDM J -t 1 ' 

1 1 • » • ‘ 

J , M 1 • 

= ■10.000 

■ ■ ' 

GAMMA = 1. 

400 



0SLTA2 

r»2/f’0 

] M? 

'-17 

.M4 

Dt;LTA3 

nELTA4 

Sir>KA2 

S rG-'A.13 

SI0MA4 

P2/P1 

Pi*/P2 

04 / Pi 

(P4/P1 1 

A 

19,0 ?. 0 

0.000 

4.? ‘19 

2.753 

0.3P5 

1 9. 52 0 

0.0 

2.5.240 

30.951 

90.000 

21.047 

5.535 

116.500 

68.263 


20,700 

0.001 

4 .07-. 

2 » 6 ft 4 

0.365 

1 3 , 97 9 

1.721 

26 . 664 

3'i .042 

69.635 

23.329 

4.004. 

116,495 

68,269 


2’ .^00 

O.OOi 

3 .9 34 

2.640 

0. 7R6 

1 8.66 5 

2.775 

27.516 

??. .132 

89.419 

24.736 

4.709 

.116.488 

66.279 


22,1 00 

0. 001 

^ r> - 4 

7. “^9 6 

0.7 P.o 

18 . 753 

3.747 

26,374 

31.249 

69.203 

26.182 

4.449 

116.477 

68.294 


??.aoo 

0.001 

3 . 7 1 P 

2.550 

C .790 

18.044 

•4,756 

29.238 

8.1 ,3 91 

88 . 987 

27,667 

4.209 

1 16.464 

68.313 


.SOO 

0-001 

3 . 60o 

2.503 

0.?OQ 

17.735 

5.765 

30.1 06 

31 .5 56 

88.770 

29.191 

3.989 

1.16.446 

68.337 


24.^00 

O.OOi 

3.4°>i 

7.45 6 

0.302 

17.426 

6 . ? 74. 

30.984 

31.751 

■68.55,1 

30.752 

3.786 

116,425 

68,365 


£'■‘.900 

O.OOl 

3 . 3 “ 3 

,2.40 5 

0.793 

17. U6 

7.78'+ 

31 . 567 

31.968 

8R. 331 

32.351 

3.5 93 

i 16.401 

6P.39P 


2‘"..d00 

O.OOl 

3 . 2 1 i 

2.760 

0.304 

1 6.604 

8,796 

32. 755 

32. 2 U 

88.108 

33. 186 

3.424 

116.373 

68,437 


2 6 . "> 0 0 

0.00'! 

3.3. '•'1 ? 

?. 712 

0.3 96 

1 0.459 

9.811 

33.651 

3 ?. 479 

87.683 

35.653 

3,263 

116,341 

6 8 . A 8 .1 


£7,000 

O.OOi 

3 .O'-- 7 

?. 26 3 

0.3 95 

16.170 

10.8 1*0 

34.554 

32.773 

87.654 

37. =AS 

3.113 

1 16.304 

63, 531 


2 7.7 00 

0.001 

3 .004 

2..215 

0.400 

15. .64 7 

]. 1 . 6 5 3 

35, 465 

33.095 

87,42.1 

39.1 OB 

2.073 

116.264 

68.587 


2!^ .^00 

0.001 

.? , ? 1 4 

?. lAA 

0.40 5 

15.519 

1 2.8BI 

36,384 

33 .445 

87,184 

40.886 

2.84.3 

116.218 

68.649 


2 •■’.I 00 

0.001 

2 • <1 ? 

2. il 7 

0.40 5 

.1 5.16 5 

13.915 

37.311 

33.824 

66.94? 

42.608 

?. 7 21 

116. 1 68 

68,71 9 


2 9. J.UO 

O.OOi 

2.741 

?. C69 

0.405 

14.844 

1 4.0 56 

,33.24,6 

34.234 

86,695 

44 , 545 

2.607 

116.112 

68.796 


''O.i^OO 

O.OCl 

? .337 

.2.0?0 

0.4 U 

14 . 49 5 

16.005 

39 . 194 

34.A77 

66.441 

46. 4 ? 5 

?,500 

1 16.050 

68.8R2_ 


31 .?00 

O.OOl 

2 .374 

i. 97? 

0.415 

14. 136 

17.062 

40.151 

35,156 

86. 160 

4R.340 

2.309 

1.15,0 32 

68.076 


31 -900 

0.001 

2 . 49 A 

1.923 

0.419 

13.771 

1 B.129 

41 .119 

35.673 

85..910„ 

50.289 

-.? + .3n5.. 

-1 1,5.007. 

-6 9,08I_ 


2 ? . M)0 

0.001 

2 .41 n 

1.575 

0.A7? 

13.394 

1 9.206 

42. 100 

36. 231 

85.632 

52.272 

2.216 

115.823 

69.197 


^3.300 

0.001 

2.3+3 

1 .5.2 7 

0-425 

13.Q04 

-20.396_ 

.4 '3 , 0 . 95 . 

.36,. 3 34 . 

„9 5.3:43. 

.54,291„2.i 32-JL15.731. 

69.826 

a 

2./-.00 0 

O', 001 

2.26 6 

1.779 

0.533 

12.601 

21.399 

44. 105 

37.437 

85.043 

56.344 

2.05? 

115.629 

69.469 

1 

700 

0.001 

? . 1 9 ? 

1.731 

0.430 

12.183 

2 ;i .517 

^.5..l'3L 

_3fi-..193..J3 4.+2Z9 

Jj8-4.24. 

-1+5.77 

-.1.15.,. il5. 

-6 9. 62a_ 


3 3 , 4 0 \) 

0.001 

7.31? 

1.453 

0 • ^ ‘I’H 

1 1.7-49 

2 3.551. 

46.177 

38. 9f.^, 

84,400 

60-562 

1.9 0 5 

115.389 

69.806 

1 

7.-, , 100 

0 ,00 

2.04 7 

1 . A 1 5 

_a+.6_?.L_.Ll . 2..?.,5. 

..?.i.8 05_ 

. 47 . 2 . 43 .. 

_39.,9 07. 

.84,053. 

62,729.. 

_1.327._ 

-115.247. 

..70-006- 


3 3, '^00 

0.002 

1 ,9 7A 

1 . 5 R 7 

0,459 

1.0. 32 0 

25.980 

48 , 334 

40.729 

63.684 

64,039 

1 .772 

1 1 5.088 

70.232 

a 

37. 300 

0.002 1.903 

1. 5 

. 0 ..46.1. 

_.L0...31 9„ 

.2.7-lSl._ 

_49,452_ 

.41. .746. 

_83.291. 

.6.7,195. 

-I.. 710. 

-114. .9 03. 

-70.489.- 


33.200 

0.00 2 

1 .2? 4 

1 .491 

0.476 

9. 790 

23.410 

50.603 

42.875 

62.869 

69,502 

1 ,650 

114.702 

70. 785 

a 

' 3.000 

0.002 

1 .7^ 7 

1 .<+4? 

0.457 

9 , 2 ? 6 

_Z.9t6 74 

5.L. 79 ? 44.1.40 

8 2,. 41 6 . 21 .. 5 9 '’•_j,i.4. 4 a 4 . 

_7±. iZS 

a 

? 3 , 0 0 

0.002 

1 .6 '-'2 

1,30? 

0.499 

8,620 

70.9B0 

53,027 

45.5 71. 

81.905 

74.200 

1.537 

U4.137 

71.532 


'^0. 3 00 

0.002 


..’4.2_ 

0.514 2,OhZ_ 

32.338 

,54. 320. 

47.214 

.81.342- 

,.76.012, 

..1.482 

1 1’.85a_72.qi 6 


41 .0 00 

0 .00? 

1 . 54 4 

1.291 

0.531 

7.236 

33.762 

55.687 

49.1 3,9 

P0.700 

79,426 

1.428 

U3.453 

72.614 


4 j . 7 () 0 

0.00 2 


1.7 '3 7 _0.5 52_ 

A. 42.6. 

_35.2?4._ 

. 1 : 7 . 153 

.51.458. 

.J. 9 . 94 .I. 

-a2.17 7„ 

1-.374. 

.212.9 45. 

.73.376 


47.400 

0.002 

1.336 

1.12 2 
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Table 14 Actual Results of Iteration from Approximate Method for Estimating Mach Stem Height, 

>' !:i'1.9‘2 > 'P' -7 p" ■= 3.67 

1 ’ ol a 


St 

Y 

o 

X 

e 

^2 ■ 

' ^3a 

4a 


"4a 

W(Y ) 
o 


0.99900 

0.00105 

0.00076 

0.00078 

0.89010 

0.00009 

0.08237 

0.04175 




0.00100 

1.04689 

0.75505 

0.76110 

0.00089 

0.08990 

0.00008 

-0.01086 



20 

0.20700 

0.83101 

0.59936 

0.60549 

0 . 18445 

0.7128 

0.01707 

0.00059 




0.19645 

0.84207 

0.60733 

0.61346 

0.17505 

0.07223 

0.01620 

0.13695 

X 

10"^ 


0.19642 

0.84210 

0.60735 

0.61348 

0.17502 

0.07223 

0.01620 

0.70090 

X 

IQ-Il 


0,99900 

0.00105 

0.00076 

0.00078 

0. 89010 

0.00009 

0.08237 

0.04175 




0.00100 

1.04689 

0.75505 

0.76138 

0.00089 

0.08993 

0.00008 

-0.01110 



40 

0.21059 

0.82725 

0.59664 

0.60300 

0.18765 

0.07098 

0.01736 

0.00060 




0.19984 

0.83852 

0.60477 

0.61113 

0.17807 

0,07195 

0.01648 

0.13036 

X 

10"^ 


0.19982 

0.83854 

0.6,0479 

0.61114 

0.17805 

0.07195 

0.01648 

0.57496 

X 

o'l^ 


0.99900 

0,00105 

0.00076 

0.00078 

0.89010 

0.00009 

0.08237 

0.04175 




0.00100 

1.04689 

0.75505 

0.76155 

0.00089 

0.08995 

, 0.00008 

-0.01124 



100 

0.21271 

0,82503 

0.59504 

0.60153 

0,18954 

0.07080 

O'. 017 54 

0.00061 




0.20185 

0.83641 

0.60325 

0.60974 

0.17986 

0.07178 

0.01664 

0.12579 

X 

lo'^ 


0.20183 

0.83643 

0.60327 

0.60976 

0.17984 

0.07179 , 

0.01664 

0.50140 

X 

10"^^ 


Experimental result = 0.20 by A. Ferri [43]. 
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Table 15 Actual Results of Iteration from Approximate Method for Estimating Mach Stem Height, 

M, =■ 1.92 , P ,/P =3.54 

1 ’ ol a 


^st 

Y 

0 

X 

e 

«2 

. 

«3a 

«4a 

r»3a 

Jo ^^3a 

jo <•••> 

W(Y^) 

0 


0.99900 

0.00102 

0.00075 

0.00077 

0.88080 

, 0.00009 

0.08112 

0.04055 • 


0.00100 

1.01372 

0.74682 

0.75214 

0'.00088 

0.08874 

0.00008 

-0.01276 

20 

0.23384 

0.77136 

0.56827 

0.57367 

0.21152 

'0,06744 

0.01947 

0.00082 


0.22884 

0.78253 

0.57650 

0.58190 

0.20182 

0.06842 

. 0.01858 

0.20009 X lO"® 


0.22880 

0.78256 

0.57652 

0.58192 

0.20178 

0.06842 

. 0.01858 

0.19920 X 10"^° 


0.99900 

0.00102 

0.00075 

0.00077 

0.88080 

0.00009 

0.08112 

0.04054 


O.OOlOO 

1.01372 

0.74682 

0.75239 

0.00088 

0.08877 

0.00008 

-0.01297 

40 

0.24282 

0.76834 

0.56605 

0.57163 

0.21415 

0.06719 

0.01971 

0.00063 


0.23164 

0.77968 

0.57440 

0.57999 

0.20429 

0.06819 

0.01881 

0.19330 X 10"^' • 









-10 


0.23161 

0;77S72 

0.57443 

0.58001 

0.20426 

0.06819 

0.01880 

0.17200 X 10 


0.99900 

0,00102 

0.00075 

0.00077 

0.88080 

0.00009 

0.08112 

0.04054 


0.00100 

1.01372 

0.74682 

0.75254 

0.00088 

0.08879 

0,00008 

-0.01309 

100 

0.24458 

0.76655 

0.56473 

0.57042 

0.21570 

0.06705 

0.01986 

0.00064 


0.23330 

0,77800 

0.57316 

0.57886 

0.20576 

0.06806 

0.01894 

0.18865 X 10~^ 


0.23327 

0.77803 

0.57318 

0.57888 

0.20573 

0.06806 

0.01894 

0.15590 X 10"^'^ 


Experimental result = 0.40 by A. Ferri [43]. 




Table 16 The Effect of Chatiging the Numb'er of, Divisions on the Initial Characteristic Line 

. on the Calculatecl Mach Disc Result 


P /P 
e a 


N 

^sk 

^sk 

^sk 

’^sk max . 

^sk 

max 

Xf at 
max 

max 


md 

20 

1.8127 

1.1523 

2.8754 

1,2241 

‘ 3.8859 

1.8961- 

5.-4369 

0.7478 

30 

1.7652 

'1.1460 

2,9067 

1.2208 

3.8960 

1.8962 

5.4216 

0.7421/ 

40 

1.6729 

1.1345 

2.9082 

1.2206 

3.8958 

1.8956 

5.4190 

0.7333 

30 

1.9097 

0.9515 

2.1042 

0.9542 

3.4187 

1.6359 

4.8779 

0..3406 . 

40 

1.8641 

0.9503 

2.1173 

0.9539 

, 3.4185 

1,6359 

4.8847 

0.3318 

25 

1.6359 

■1.36,57 

3.7906 

1.6762 

4.6869 

,2.3416 

6.1909 

1.3478 

40 

1.4843 

1-.3273. 

3.9601 

1.6695 

' 4.7881 

'2.3429 

6.1838 

1.3444 

25 

2.1300 

1.3192 

4.6932, ■ 

1.5623 

6.0133 

2.2995 

8.2792 

0.9055 

40 

1.9856 

1.2930 

4-. 552 8 

1.5573 

5.9322 

2.2998 

8.3509 

0.8604 

25 

2.2450 

1.4686 

6.9194 

2,0842 

8.4816 ' 

2.8652 

12.4430 

1.1712 

40 

2.0939 

1,4322 

7.0029 

2.0780 

8.5051 

2.8736 

12.5155 

1.1558 

30 

2.4854 

1.3249 ' 

5.6870 

1.6130 

7.6451 

2.3909 

10.8680 

0.6680 

40 

2.3525 

1.3038 

5.7954 

1.6099, 

7.5325 ' 

2.3894 

11.0602 

0.5718 

30 

2.6189 

1.3929 

7.8606 . 

1.8927 

9.8373 

2.7064 

14.9062 

0.5927 

40 

2.6163 

1.3911 

7.7557 

1.8941 

■ 9.8856 

2.7045 

14.9046 

0,5707 

20 

2.7939 

1. 5110 

8. ‘6360 • • 

2.1330 

10.3414 

2.9484 

16.207 

0.7031 

. 30 

2.5454 

1.4592 

8.3920 

2.1206 

10.3385 

2.9421 

15.8100 

,0.8307 


Oi 

o 
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APPENDIX A 

FUNCTIONS FOR EQS. (24) and (25) OF TWO-DIIffiNSIONAL EXTERNAL 


See Fig. 13a and References ‘[1? ,25 ] . 
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2 ds o so 
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(A. 2) 

(A. 3) 
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(A. 7) 


(A. 8) 
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and 

K -*■ 


K, 


" ["t (^ ■ 2 ' S ' ' ) ~ 

L ^ M cos K ' 


cos 0 - cos K sin K sin 0 (A. 9) 


•^3 “ 2 (^ ■ 2 2 )“ SI 

' M cos K ' 

QO 


. 2 
sin k: 


sin 6 + cos K sin K cos 6 (A. 10) 




J 

f2 + (y - 1) cos^ K~] 

' CO 

y 

y(y - 1) tan k 

Y + 1 

(y + 1) cos^ K 

^ CO ^ 


2+(y-l)M^ cos^ K 

U * CC J 




cos^ K - ij 
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cos 1C - — ^ 

M M 


)■ 


cos 6 


2 

Y + 
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(A. 13) 


r = u 
ne - “ 




COS^ < - 


) 


(-sin 0) 


2 ^ 2 

1 \1 

cos 0 

- tan icicos: k - 

2 / 

J X X \ 


J 


ly(y-l) - [l/(y-D] 

p = (1 - W^) (|) 


P = 4* 


-I 


, 2 2 
(i) = / V + V 

s n 


(A.14) 

(A. 15) 
(A. 16) 
(A. 17) 


for attached shock wave, K -> ic as £ -^ 0; for dexached shock wave 

o . 


0 as E -> E . 

o 


*In Reference [17] for axisymmetric ’ case , J = 2. Hovrever, in the present 
derivation j is twice as large, that is , J = 4 for -both axisymnietric 
and two-dimensional cases, which agrees with the results of [25]. 
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FUNCTIONS FOR EQS. (35) and (36) ‘OF TWO-DIMENSIONAL INTERNAL FLOW 


See Fig. 13b and Reference' [21] . 
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where prime (') indlcetes the derivative taken with respect to the 
horizontal coordinate x. 



APPENDIX C 


CHARACTERISTIC EQUATIONS FOR NUMERICAL CALCULATION OF A 
STEADY SUPERSONIC FLOW FIELD 

The character»is1:ic equations fox> a ti-yo-dimensional or an axisym- 
metric flow are written as [31] : 


dr 

dx 


tan (6 + a) 


I, II 


(C.13 


. 1 dM^ 


d 0 


= + tan a + 


V tan a tan 6 dr 


I, II 


tan 0 + tan a 


rd0 


sin 




where v =■ 


Y(Y - 1) d0 

-0 for two-dimensional flov? 

'1 for axisymmetric flow- 
■0 for irrotational flow • 


I, II 




CC.2] 


•1 for rotational flow 
Furthermore , for an adiabatic process of a perfect gas , the change 
in entropy can be related to the variation of stagnation pressure 
through the second law of thei^.odynamics . 


r -Cy - 1) In 5 — 

L^olJ 


(C.3) 


C.l , PRANDTL-MEYER CORNER EXPANSION 

At ‘the comer point where the expansion from state 1 to state 2 is 
two-dimensional and isentropic, that is , V = 0, 52 = 0, the direct inte- 
gration of Eq. (C.2-) will give 
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6 ^ [ = + CwCM§) - u(H*Oj 

[I, II 


(C.4) 


vrhere 


tt)( M* ) 




T I 
Y - 1 



1 





Y 


- 1 


Y + 1 




(C.5) 


is the Prandtl-Meyer function. 

Equation (C.4) is used to establish the flow properties of the ex- 
pansion fan at the nozzle exit comer or. at the location where shock 
intersects X'rith the free j'et boundary. 


C.2 FIELD POINT PROCEDURE 


As shown in Fig, C.l, points 1 and 3 are joined by a characteris- 
tic line of family I and points 2 and 3 are joined by a characteristic 
line of family II. - Points 1 and 3 are knovm, and point 3 is to be de- 
termined. Equations (C.l) and (C-.2)j.when written, in finite difference 
-forms , give 


^3 ^ 

^3 " ^1 


tan ( ) 


(C. 6 ) 


(023 ^ « 23 ^ 


(C.7) 


A ^^3 - 

tTa 63 ~ 6 - 

•^l3 ^ i 


= -tan a 


V tan tan 6^3 


^3 " ^1 


13 * — 

tan 6^.3 + tan a 
Xo 


(e - e ) 


13 13 3 


n . 2 — ■ 

n sin a 

^ Vc ' e3 ' ep ^ 


o3 


ol 


(C. 8 ) 
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- - ^ ^3 ^23 "^3 " "^2 

tan a +. ^ 2!^ 

tan 023 + tan a23 r2s ^®3 “ ®2'^ 


. 2 ~ 

0 Sin «22 



(C.9) 


where Eq. CC.3) has-been applied. The bar superscript indicates that 
an average value is to be taken . 

In the small characteristic grid shovm in Fig. C.l, if the stream- 
line passing through point 3 intersects line 12 at point s, the co- 
ordinates of point s can be found from the geometric relation of the 
intersection of two straight lines, 12 and s3. Hence, 
r - r + X tan 0 - x tan B 

X = 2 -3 ? ^ ^ .■ 

^ - tan 6 - tan 3 

S ' * 


_ . r^ = t '2 - -.x^) tan 3 (C.ll) 

The flow angle and stagnation pressinpe at point s can , then , be approxi- 
mated by using linear interpolation. 




(0 




(C.12) 


p =p =p -h(P — •— = — 

os o3 o2 ol o2^ r^^ - r 2 


where P = P „ - because point 3 and point s are on the same streamline 
os o3 ’ 

■and have the same entropy. 

The iterative solution of the simultaneous equations, Eqs . (C.6), 
(C,7), (C.S), and (C.9), with the help of Eqs. (C.IO), (C.ll), (C.12), 
and (C.13), ‘.vill furnish the unknovm flovr variables at point 3. 

' To begin with an iterative solution, the property at point 1 is 
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usually used as the average value between points 1 and 3, and that at 
point 2 is used between points 2 and 3. However, Eqs. (C. 8 ) or (C.9) 
will blow up when points 1 or 2 are on the axis of symmetry. This 
difficulty can be avoided through the following modification: Refer- 

ring to Fig. C.-2, when point 2 is on the axis of symmetry, and intro- 
ducing • 


lim 


tan tan 


tan 0 


r^ 0 tan + tan ^ 23^2 3 

02-0 


lim 


23 


^2 " ° ®23 


Vtan a^g / 


•23 


because lim 

. ej » 0 


®23 “ ®23 


23 

’23 


(C.14) 


Equation ( C . 8 ) becomes 


Hi 




M* 

23 


= 63(1 +.v).tan «23 ^ 


Y \>’o2 


(C.15) 


which is used to replace Eq.' (C.9) for the special case when point 2 
Is on the axis of symmetry. 


e.3 CONSTANT PRESSURE JET BOUNDARY POINT PROCEDURE 

1 As shown in Fig. C.3, points 1 and 3 are on the same free jet 
boundary without intersecting with shock discontinuity, and points 2 
and 3 are joined by a characteristic line of family II. Points 1 and 

2 are known, and point 3 is to be calculated. 

The slope of the jex boundary between points 1 and 3 is given by 



tan 0 


13 


(c;l6) 
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On the constant pressure jet boundary 



and’, since free jet boundary is also a streamline, 

P , = P , 
o3 ol 


(C.17) 


(C.18) 


Therefore, Eqs. (C.7), (C.9), and CC.15), plus conditions (Eqs. (C.17) 

f 

and (C.18)) can be solved for x , r , and 6 . 

O O o 


Q.4 AXIS POINT PROCEDURE 

■ If point 3 is at the intersection of a characteristic line and 
the axis of symnetry as shown in Fig. C.3 and points 1 and 3 are 
joined by characteristic curve of family I. Points 1 and 2 are known, 
and point 3 is to be determined. 

Since points 2 and 3 are on the same axis streamline, 

■ '■<>3 =-^02 
Also, on the axis of symmetry 

= 0, - Tg = 0 (C.20). 

Therefore, Eqs. (C.6) and (C.8) can be used to solve for tvfo unknowns, 
K| and Xg. 


C.5 DOWNSTREAM POINT OF A SHOCK WAVE 

Downstream of a shock wave , the points 1 and 2 are known ; point 3 
just behind the shock v?ave is to be determined. Referring to Fig. C.5, 
points 2 and 3 are joined by characteristic line of family II. In front 
of the. shock wave , the flow is nonuniform and the flow properties are 


known. 
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The slope of 'the shock wave bet^^een points 1 and 3 is given by 

(C.21) 


^3 ■ ^1 ^ -jr 


' 3 - 


where 


^ ^ ®i»> * 'PsI + 

^13 - 2 


(C.22) 


Deflection' angle 5^ and stagnation pressure are readily available 

from oblique shock wave equations whenever the shock wave angle is 

assumed knovm during the iteration process (Note: 0^ = 

Equations (C-21), (C.7), and (C.9), together v/ith the oblique 

shock wave relations, are sufficient to determine the unknowns, x^, r^, 

M*, 8„, and P „ at point 3. 

3 3 oo 
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Figure C.3 Jet Boundary Point 



Figure C.2 Field Point with 
Point 2 on Axis 
of S 3 nn. 



Figure C.4 Axis Point 



Figure C,5 Downstrean Point of a Shock VJave 
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